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ENTRAINMENT  DESCRIPTIONS  FOR  MATHEMATICAL  MODELING 


OF  PUMPED- STORAGE  INFLOWS  IN  RESERVOIRS 


PART  I :  INTRODUCTION 


Background 

1.  Pumped-storage  (PS)  hydropower  is  a  well  established  means  of 
storing  large  amounts  of  energy  for  subsequent  power  generation.  Pumped 
storage  projects  pump  water  to  an  upper  reservoir  during  periods  of  low 
energy  demand  when  power  can  be  bought  at  a  reduced  price.  The  water 
stored  in  the  upper  reservoir  is  later  released  during  periods  of  peak 
power  demand  to  generate  high-value  power. 

2.  Prior  to  construction  of  a  PS  project,  the  impact  of  the  proj¬ 
ect  upon  existing  water  quality  should  be  assessed.  The  PS  inflows  and 
outflows  cause  mixing  that  can  significantly  alter  thermal  stratifica¬ 
tion  and  the  water  quality  parameters  within  the  reservoir.  Conse¬ 
quently,  distributions  within  the  reservoir  affect  the  temperature  and 
quality  of  downstream  releases.  The  effect  of  PS  inflows  is  well  demon¬ 
strated  by  the  temperatures  observed  at  Allegheny  Reservoir,  Pa.,  for 
pre-  and  post-PS  hydropower  (Dortch  1981). 

3.  Comprehensive  modeling  of  PS  projects  necessitates  simulation 
of  a  complex  dual  reservoir  and/or  riverine  system.  The  quality  of  in 
situ  and  release  water  from  conventional  reservoirs  is  dependent  upon 
meteorology,  hydrology,  internal  mixing  processes,  morphometry,  regu¬ 
lating  structures,  and  operational  methodology.  The  PS  system  is  more 
complex.  One  of  the  reservoirs  of  a  PS  system  is  usually  much  larger 
than  the  other  and  normally  goes  through  a  thermal  stratification  cycle, 
whereas  the  smaller  reservoir  usually  does  not  stratify  vertically. 
Relatively  large  amounts  of  water  are  passed  between  the  two  reservoirs 
during  generation  and  pumpback,  thus  causing  mixing  and  redistribution 
of  thermal  and  water  quality  characteristics. 

4.  When  the  PS  flow  enters  the  larger  stratified  reservoir 


(Figure  1),  a  jet  is  often  produced  that  entrains  ambient  water.  The 
entrained  ambient  water  combines  with  the  entering  flow  and  causes  an 
increase  in  the  volume  flux  of  the  jet.  Density  stratification  of  the 
ambient  pool  eventually  causes  the  jet  to  collapse  and  change  into  a 
tranquil  density  current  without  additional  entrainment.  This  entrain¬ 
ment  process  can  significantly  alter  the  temperature  (Figure  2)  and 
water  quality  parameters  in  the  water  column.  Previous  studies  (i.e. 
Dortch  1981)  have  demonstrated  that  the  entrainment  description  can  be 
one  of  the  most  important  inputs  to  a  PS  reservoir  simulation  model. 
Therefore,  adequate  entrainment  descriptions  are  necessary  for  realistic 
models  of  PS  projects. 


Figure  1.  Schematic  of  pumped-storage  inflow  jet 
and  entrainment 


Past  Experience 

5.  A  hybrid  technique  has  been  used  successfully  at  the  US  Army 
Engineer  Waterways  Experiment  Station  (WES)  to  model  PS  hydropower  proj¬ 
ects  on  a  variety  of  existing  and  proposed  PS  projects  including 
Richard  B.  Russell  Reservoir  (Smith  et  al.  1981,  Fontane  and  Bohan 
1974),  Marysville  Lake  (Fontane  et  al.  1977,  Dortch  1978),  Dickey- 
Lincoln  School  Lakes  (Dortch  et  al.  1976),  and  Kinzua  Reservoir  (Dortch 
1981).  This  approach  employs  a  one-dimensional  (1-D)  (vertical)  reser¬ 
voir  simulation  model  that  includes  algorithms  for  the  pumped-storage 
mixing  process.  The  algorithms  were  developed  on  a  site-specific  basis 
from  observed  PS  inflow  characteristics  produced  in  physical  models. 
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6.  The  primary  information  obtained  from  physical  models  regarding 
PS  inflows  was  the  entrainment  description.  The  entrainment  description 
prescribes  the  total  entrainment  rate  E  (as  a  multiple  of  inflow  rate) 
and  the  vertical  distribution  of  entrained  flow  Q&(z).  These  param¬ 
eters  were  obtained  from  the  physical  model  for  a  typical  set  of  flow 
and  stratification  conditions,  and  the  values  were  used  throughout  the 
numerical  simulations.  This  procedure  did  yield  fairly  accurate  tem¬ 
perature  simulations,  as  demonstrated  by  Figure  3  from  the  Allegheny 
Reservoir  Study  (Dortch  1981). 

7.  However,  the  entrainment  parameters  are  actually  not  constant 
over  a  simulation  period.  These  parameters  are  dependent  upon  the 
source  geometry,  buoyancy,  momentum,  and  volume  flux  and  the  ambient 
stratification.  Therefore,  entrainment  can  vary  from  day  to  day  at  a 
particular  project  for  different  stratification  and  flow  conditions. 

The  temperature  distribution  is  certainly  influenced  by  the  entrainment 
description,  but  the  simulated  distribution  has  been  found  to  be  mildly 
sensitive  to  variations  in  the  entrainment  parameters,  especially  E  . 

For  example,  for  the  Allegheny  Reservoir  Study,  the  temperature  results 
were  found  to  be  relatively  insensitive  to  variations  in  E  of  up  to 
about  50  percent.  However,  it  was  important  to  use  the  correct  propor¬ 
tions  of  entrainment  from  the  epilimnion  and  hypolimnion.  Other  param¬ 
eters  that  are  less  continuous  in  the  water  column  than  temperature  may 
be  more  sensitive  to  the  prescribed  amount  of  entrainment  and  dilution. 


TEMPERATURE,  *C  TEMPERATURE,  "C  TEMPERATURE.  ”C 

11  JUL  1977  8  AUG  1977  12  SEP  1977 

Figure  3.  Comparison  of  simulated  and  observed  temperature 
profiles  at  Allegheny  Reservoir  (Note:  1  ft  =  0.3048  m) 


Thus,  for  reservoir  water  quality  simulations  of  PS  projects,  it  seems 
desirable,  if  not  necessary,  to  have  entrainment  algorithms  that  are 
functionally  dependent  on  time-varying  flow  and  stratification 
conditions . 

Purpose  and  Scope  of  Study 

8.  The  purpose  of  this  work  was  to  develop  general-purpose  algo¬ 
rithms  to  predict  the  entrainment  parameters  for  varying  types  of  PS  in¬ 
flows  with  varying  conditions.  The  algorithms  had  to  fit  into  the  frame¬ 
work  of  1-D  reservoir  water  quality  models,  such  as  CE-QUAL-R1  (Environ¬ 
mental  Laboratory  1982).  It  should  be  realized  that  general-purpose 
entrainment  algorithms  cannot  be  used  for  every  project  without  some  ad¬ 
justments.  For  example,  local  topography  and  unusual  source  geometries 
can  influence  entrainment  and  mixing.  In  such  instances,  the  hybrid 
(physical-numerical  models)  approach  discussed  earlier  is  recommended 
either  to  select  and/or  calibrate  the  most  appropriate  entrainment  de¬ 
scriptions  or,  if  necessary,  to  develop  alternate  descriptions.  In 
either  case,  a  general-purpose  algorithm  is  more  expedient  for  defining 
the  appropriate  description. 

9.  Following  a  literature  study  by  Roberts  (1981),  it  was  con¬ 
cluded  that  existing  jet  models  could  be  used  to  describe  most  of  the 
highly  buoyant  types  of  PS  inflows  that  might  be  encountered.  However, 
little  or  no  information  was  available  for  the  prevailing  cases  of  weak 
buoyancy.  Pumped-storage  inflow  often  takes  the  form  of  a  jet  of  fairly 
high  momentum  but  small  buoyancy  discharging  into  a  density-stratified 
ambient.  There  are  few  experimental  data  available  on  this  class  of  jet 
flows.  As  part  of  this  study,  a  comprehensive  series  of  experiments  on 
jet  configurations  typical  of  PS  inflows  was  performed.  These  include 
horizontal  and  inclined  jets  of  neutral  or  small  positive  or  negative 
buoyancy  discharging  into  stratified  ambients.  These  experiments  and 
the  resulting  algorithm  for  weakly  buoyant  jets  are  described  herein. 

10.  Along  with  the  weakly  buoyant  jet  algorithm,  a  jet  model  for 
strongly  buoyant  inflows  was  selected  and  incorporated  into  the  1-D  res¬ 
ervoir  simulation  code.  In  addition  to  the  discussions  of  the  weakly 
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PAST  II:  MODELING  FRAMEWORK 


11.  Mathematical  models  have  been  used  in  the  past  (Dortch  et  al. 
1976,  Dortch  1981,  Holland  et  al .  1982)  to  simulate  thermal  stratifica¬ 
tion  and  release  temperatures  in  PS  (as  well  as  conventional)  reservoirs. 
Previous  studies  had  indicated  that  for  PS  reservoirs,  the  PS  inflow  was 
a  dominant  mixing  process.  Therefore,  these  models,  which  included  nat¬ 
ural  stratifying  and  mixing  mechanisms  such  as  surface  heat  inputs,  in¬ 
flow,  outflow,  wind  mixing,  and  internal  diffusion,  were  modified  to 
include  mixing  mechanisms  for  PS  operations. 

12.  A  1-D  (vertical)  reservoir  thermal  simulation  code  (WESTEX) , 
which  was  developed  by  WES,  was  modified  to  account  for  PS  inflows  and 
entrainment.  The  code  also  allows  simulation  of  an  afterbay,  which  of¬ 
ten  is  a  part  of  a  PS  project.  The  larger,  stratified  reservoir  of  the 
system  is  conceptualized  as  a  number  of  homogeneous  horizontal  layers 
stacked  vertically,  and  the  heat  sources  and  sinks  to  a  general  layer 
are  represented  as  shown  in  Figure  4.  If  a  dual-reservoir  system  is 
used,  then  the  other  reservoir  must  be  modeled,  too.  In  all  previous 
studies,  the  smaller  reservoir  was  not  expected  to  stratify;  thus,  fully 
mixed  conditions  were  used  to  budget  heat  and  water. 

13.  The  solution  for  the  temperature  history  of  a  general  layer 
is  obtained  through  conservation  of  mass  and  energy  and  an  equation  of 
state  for  density  as  a  function  of  temperature.  The  energy  equation  is 
stated  as: 
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An  equation  similar  to  Equation  1  would  exist  for  each  water  quality 
constituent. 


14.  The  pumped-storage  inflow  Q 


Q^(z)  from  all  layers  and  the  PS  initial  jet  inflow  rate 


is  the  sum  of  entrained  flow 

The  PS 

inflow  temperature  0^  is  a  flow-weighted  average  of  the  inflow  jet  at 
its  temperature  and  all  the  entrainment  flows  at  their  respective  layer 
temperatures.  The  temperature  of  the  inflow  jet  is  obtained  from  the 
heat  budget  of  the  afterbay  or  the  smaller  reservoir  of  dual-reservoir 
systems . 


15.  The  entrained  flow  is  treated  as  withdrawal  of  specified 
amounts  from  specified  layers.  It  is  important  to  realize  that  entrain¬ 
ment  withdrawal  results  in  vertical  flows  downward.  The  PS  inflow, 

Qp  ,  is  handled  in  the  same  manner  as  the  upstream  inflow;  it  is  placed 
into  the  layer  of  neutral,  or  closest  to  neutral,  buoyancy.  The  inflows 
cause  vertical  upward  flow.  The  entrainment  distribution  and  rate, 
which  are  affected  by  source  geometry,  buoyancy,  volume,  and  momentum 
fluxes  and  ambient  stratification,  must  be  specified.  The  total  en¬ 
trained  flow  rate,  ,  can  be  expressed  in  terms  of  the  jet  volume 

flux  through  a  total  entrainment  coefficient,  E  , 

QE  =  EQj  =  2  Qe(z)  (2) 

all  layers 


resulting  in 

Qp  *  O  +  E)Q 


(3) 


If  the  shape  of  the  entrainment  distribution  is  known  or  specified,  the 
entrained  flow  to  be  withdrawn  from  each  layer  (Qg(z))  can  be  determined. 

16.  A  daily  time  step  for  the  simulations  has  been  found  to  be  ade¬ 
quate.  However,  the  PS  inflow  rates  must  be  used  to  obtain  entrainment 
flows  that  are  functionally  dependent  on  momentum  and  volume  fluxes. 
Additionally,  outflow  rates  are  required  to  predict  the  withdrawal  zone. 
Since  PS  hydropower  outflow  and  inflow  rates  can  vary  significantly  dur¬ 
ing  a  daily  time  step,  actual  release  and  pumpback  rates  should  be  used 
with  the  durations  of  these  flows  to  account  for  changes  in  daily  water 
budget . 


PART  III:  JET-INDUCED  ENTRAINMENT 


Introduction 


17.  In  a  literature  survey  of  methods  of  predicting  entrainment  in 
the  types  of  jets  expected  in  PS  reservoirs,  Roberts  (1981)  identified 

16  possible  jet  configurations  that  could  arise.  Of  these,  some  had 
been  well  studied  and  were  well  understood,  and  some  had  not  been  stud¬ 
ied  at  all.  Subsequently,  it  was  decided  to  investigate  further  the 
types  of  jets  that  would  be  expected  to  occur  most  often  in  PS  reser¬ 
voirs.  These  are  horizontal  or  inclined  discharges  of  positive,  nega¬ 
tive,  or  neutral  buoyancy  into  an  arbitrarily  stratified  environment. 

Even  with  this  broad  range  of  variables  it  was  found  that  the  jets  fell 
into  two  classes:  strongly  or  weakly  buoyant.  The  former  class  has 
been  fairly  extensively  studied,  and  well-established  mathematical 
models  exist  to  predict  their  behavior.  For  this  case,  incorporation 
into  a  reservoir  thermal  model  simply  required  adaptation  of  an  existing 
code.  However,  for  weakly  buoyant  jets,  it  was  found  that  little  infor¬ 
mation  existed  of  either  a  theoretical  or  experimental  nature.  For  that 
reason  an  extensive  series  of  experiments  of  a  fundamental  nature  were 
performed  on  weakly  buoyant  jets  discharged  into  arbitrary  stratifica¬ 
tion.  The  results  were  then  reduced  to  a  form  suitable  for  inclusion 
into  reservoir  simulation  models. 

18.  The  development  of  the  algorithms  for  predicting  jet-induced 
entrainment  is  outlined  in  the  following  sections.  First,  the  general 
problem  is  discussed  and  analyzed  by  dimensional  and  length  scale  argu¬ 
ments,  and  criteria  for  separation  into  weakly  and  strongly  buoyant  jets 
are  presented.  Next,  the  algorithms  for  predicting  strongly  buoyant  and 
weakly  buoyant  jet  dynamics  are  presented.  Also,  the  results  of  the  ex¬ 
periments  on  weakly  buoyant  jets  (discussed  in  detail  in  Appendix  A)  are 
summarized,  and  the  empirical  formulas  used  to  incorporate  the  results 
into  the  reservoir  thermal  model  (WESTEX)  are  presented.  Finally,  the 
approach  is  summarized,  and  flowcharts  for  the  entrainment  subroutines 


are  given. 


19.  The  problem  under  consideration  is  shown  in  Figure  5.  A  jet 
inclined  upwards  at  an  angle  6  to  the  horizontal  discharges  into  a 
stagnant  reservoir  of  arbitrary  stratification  denoted  by  p(z)  ,  where 
p  is  water  density  as  a  function  of  vertical  position  z  .  The  jet 
entrains  ambient  water  as  it  rises  or  falls  and  may  eventually  reach  a 
height  at  which  its  density  is  equal  to  that  of  the  ambient  water.  Near 
this  height,  the  jet  collapses  vertically  and  intrudes  as  a  horizontally 
flowing  density  current.  This  study  was  particularly  concerned  with  the 
amount  of  flow  entrained  into  the  jet  and  the  vertical  distribution  of 
the  entrained  flow.  Because  of  the  large  number  of  parameters  involved 
and  the  arbitrary  shape  of  the  ambient  density  profile,  an  unreasonably 
large  number  of  dimensionless  parameters  and  experiments  would  be  needed 
to  describe  all  possible  situations.  The  problem  can  be  made  tractable, 
however,  by  making  some  simplifying  assumptions  for  the  ambient  strati¬ 
fication,  which  when  combined  with  judicious  dimensional  analysis  will 
yield  results  of  an  accuracy  consistent  with  modeling  requirements. 

20.  The  idealized  problem  is  shown  in  Figure  6.  A  jet  of  area  A 
is  discharged  at  an  angle  6  into  a  stagnant  ambient  fluid.  The  am¬ 
bient  stratification  is  characterized  by  a  well-mixed  surface  layer  sep¬ 
arated  by  a  density  jump  of  magnitude  Ap  from  a  linearly  stratified 
lower  layer.  This  stratification  is  an  idealization  of  that  occurring 
in  reservoirs.  The  idealized  problem  incorporates  several  other  prob¬ 
lems.  For  example,  the  density  interface,  if  strong  enough,  approxi¬ 
mates  a  nearby  solid  boundary.  If  the  interface  is  far  removed  from  the 
jet,  it  has  no  effect,  and  the  flow  is  the  same  as  that  induced  in  a  com 
pletely  linear  stratification.  By  linearizing  the  ambient  stratifica¬ 
tion  over  the  extent  of  jet  influence,  the  results  for  linear  stratifi¬ 
cation  can  be  used;  this  approximation  will  be  utilized  later. 

21.  The  problem  is  analyzed  by  dimensional  and  length  scale  argu¬ 
ments  in  Appendix  A  (Analysis  section).  It  is  shown  that  the  problem 
can  be  characterized  by  the  relative  magnitudes  of  five  length  scales 
that  can  be  formed  from  the  parameters  involved.  These  length  scales 


Figure  5.  Generalized  problem 


Figure  6.  Idealized  problem 


(£_,£,£  ,£A)  and  h  ,  and  their  physical  significances,  are  discussed 
Q  M  £  A 


in  Appendix  A. 

22.  Consider  first  the  case  where  h  ,  the  jet  submergence  be¬ 
low  the  interface,  is  much  larger  than  the  other  length  scales,  i.e., 


h  >>£,£,  £  ,  £  .  In  this  case,  the  jet  and  its  extent  of  in- 

Q  n  C  A 


fluence  are  far  removed  from  the  interface,  and  the  presence  of  the 
interface  should  not  affect  the  flow  field.  The  only  length  scales  are 
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now  ,  i  ,  and  £  If  £.  is  not  large  relative  to  £u  or  i  , 
y  n  £  n  € 

the  nozzle  geometry  should  exert  only  a  secondary  influence  on  the  flow 

field  and  not  affect  its  gross  characteristics  (Wright  1977b).  This 

leaves  only  £„  and  £  ,  which  characterize  the  relative  effects  of 

the  ambient  stratification  and  the  buoyancy  and  momentum  fluxes.  If 

£w  <<  £  ,  then  the  effects  of  the  source  momentum  flux  are  confined  to 

n  e 

the  region  near  the  nozzle.  The  flow  is  predominantly  buoyancy  driven, 

similar  to  that  of  a  sewage  discharge  in  a  stratified  ocean.  On  the 

other  hand,  if  £u  »  £  ,  the  buoyancy  is  small  and  the  flow  is  momen- 
n  e 

turn  dominated.  For  the  former  case  (£^  <<  £^),  mathematical  (integral) 
models  have  been  developed  and  applied  (Koh  and  Brooks  1975).  These 
models  break  down,  however,  as  £^/££  becomes  large,  and  there  is 
little  information  available  on  these  flow  types. 

23.  It  is  shown  in  Appendix  B  (see  Figure  B2)  that  the  following 
ranges  apply:  for  £M/££  <  ~°-8  ,  the  flow  is  buoyancy  dominated,  i.e. 
plumelike;  for  0.8  <  <  2  ,  both  momentum  and  buoyancy  are  impor¬ 

tant;  and  for  £^/£e  >  2  ,  the  momentum  is  dominant.  The  integral  models 
work  reasonably  well  so  long  as  momentum  does  not  dominate;  therefore, 
they  are  proposed  for  the  PS  modeling  when  £M/££  <  2  .  The  adaptation 
of  the  integral  models  to  WESTEX  is  described  in  the  following  section. 
For  £M/££  i  2  ,  the  experimental  studies  that  were  performed  as  part  of 
the  present  effort  (described  in  Appendix  A)  are  used.  The  adaptation 
of  these  experimental  results  to  WESTEX  is  described  in  paragraphs  34-38. 

Strongly  Buoyant  Jets 


Integral  buoyant  jet  model 

24.  These  jets  are  defined  as  having  £^/££  <  2  .  Several  models 
have  been  used  to  predict  these  flows;  the  well-established  model  of  Fan 
and  Brooks  (1969)  and  Koh  and  Fan  (1970)  are  used  here.  The  formulation 
of  this  model  is  described  in  Appendix  C. 

Adaptation  of  integral  buoyant 
jet  model  to  reservoir  thermal  model 


25.  The  integral  buoyant  jet  model  requires  several  modifications 


for  use  in  a  thermal  model  for  a  PS  reservoir.  These  include  choice  of 
nozzle  size,  choice  of  model  coefficients,  and  prediction  of  the  en¬ 
trained  flow  distribution.  These  adaptations  are  discussed  below.  No 
assumptions  on  linearity  of  the  density  profile  are  required,  as  the 
model  will  accommodate  arbitrarily  shaped  ambient  density  profiles.  The 
integral  buoyant  jet  model,  as  adapted  to  the  reservoir  simulation  code, 
is  referred  to  as  subroutine  SBJET  (Strongly  Buoyant  JET) . 

26.  Fan  and  Brooks'  model  assumes  a  round  nozzle  and  a  radially 
symmetric  flow  field  around  the  trajectory  of  the  jet  center  line.  For 
a  noncircular  nozzle,  the  flow  field  will  not  be  radially  circular  near 
the  discharge.  However,  as  discussed  in  Appendix  A,  for  distances  from 
the  nozzle  much  longer  than  its  characteristic  size,  the  nozzle  size 
will  not  explicitly  affect  the  flow  field.  At  this  point,  the  flow 
field  will  be  radially  symmetric.  Therefore,  we  can  approximate  the 
flow  as  being  equivalent  to  one  discharged  from  a  round  nozzle  with  the 
same  momentum  and  buoyancy  fluxes  as  the  noncircular  nozzle.  This  re¬ 
quires  that  the  equivalent  round  nozzle  has  the  same  cross-sectional 
area  A  ,  as  the  actual  nozzle,  or: 


and  d  = 

where  d  is  nozzle  diameter.  These  approximations  are  used  in  SBJET. 
In  practice,  the  approximations  will  be  better  the  closer  the  discharge 
is  to  being  round.  Thus,  a  square  nozzle  would  be  a  better  approxima¬ 
tion  to  the  equivalent  round  nozzle  than  would  be  a  long  slot. 

27.  Use  of  the  integral  model  requires  specification  of  the  en¬ 
trainment  coefficient  a  and  turbulent  Schmidt  number  \  (see  Appen¬ 
dixes  B  and  C).  It  is  well  known  that  the  entrainment  coefficient  of  a 
plume  is  greater  than  that  of  a  jet,  and  that  the  entrainment  coeffi¬ 
cient  of  a  jet  which  makes  a  transition  to  a  plume  will  vary  along  its 
trajectory.  The  entrainment  coefficient  will  increase,  being  close  to 
that  of  a  jet  near  the  nozzle  and  approaching  that  of  a  plume  far  away. 


The  values  of  the  two  coefficients  are  given  in  Appendix  B,  Equations 
B21  and  B22.  For  PS  reservoir  jets,  the  density  differences  are  small, 
and  it  is  unlikely  that  the  flow  will  ever  approach  that  of  a  plume. 

The  momentum  will  usually  dominate,  and  so  it  is  appropriate  to  use  jet 
entrainment  coefficients.  (For  further  discussion  of  this,  see  Appendix 
B.)  Therefore,  a  from  Equation  B21  is: 

a  =  0.0535  (6) 

The  term  A  is  the  same  for  jets  and  plumes;  the  value  is  given  by  Equa¬ 
tion  B23  as: 

\  =  1.2  (7) 

28.  The  most  difficult  adaptive  problem  is  the  specification  of 
the  vertical  distribution  of  entrained  volume  flux  into  the  jet,  essen¬ 
tial  for  use  in  a  reservoir  simulation  model.  This  problem  has  received 
almost  no  discussion  in  the  literature;  the  experiments  performed  as 
part  of  the  present  studies  (Appendix  A)  are  apparently  the  first  in 
which  the  entrained  flow  rate  and  its  vertical  distribution  have  been 
directly  measured.  These  experiments  do  not  cover  the  more  buoyant  jets 
for  which  the  integral  model  is  used,  however,  and  some  assumptions  are 
necessary  to  predict  the  distribution. 

29.  The  entrainment  into  the  jet  is  given  by  the  well-known  assump 
tion  (Equation  C6  and  Table  Cl)  that  the  entrained  flow  velocity  is  pro¬ 
portional  to  the  local  jet  center-line  velocity.  For  a  round  jet,  that 
is : 

^  =  2naub  (8) 

ds 

where 

3 

Q  =  jet  volume  flux,  L  /T 
s  =  distance  along  jet  trajectory,  L 
u  =  local  jet  center-line  velocity,  L/T 
b  =  local  jet  half-width,  L 
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Figure  7.  Method  for  predicting  entrained  flow  distribution 
from  integral  buoyant  jet  model 

31.  Consider  the  i*"*1  model  layer.  The  intersection  of  the  jet 
perimeter  subtends  an  angle  6  (denoted  by  ALPHA  in  SBJET)  at  the  cen 
ter.  The  perimeter  P  of  the  jet  contained  within  layer  i  is: 
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Xj  ,  X2  ’  yi  ’  an<^  ^2  anc*  t^le  vo^unie  entraine<l  are  computed  at  each 
space  step. 

33.  The  entrained  flow  distribution  for  each  layer  is  computed  for 
each  space  step  by  Equations  10-17  and  added  to  those  computed  in  the 
previous  step.  This  is  continued  until  the  terminal  rise  height  of  the 
jet  is  reached.  The  procedure  is  shown  in  the  flowchart,  Figure  8. 

Weakly  Buoyant  Jets 


Experimental  study 

34.  Weakly  buoyant  jets  are  defined  as  those  having  —  2  • 

The  extensive  experimental  program  to  investigate  these  jets  is  de¬ 
scribed  in  Appendix  A.  Adaptation  of  these  results  to  the  thermal 
model  is  discussed  below. 

Adaptation  to  reservoir  thermal  model 

35.  The  experiments  performed  on  weakly  buoyant  jets  were  of  a 
fundamental  nature,  and  some  adaptations  of  these  results  were  necessary 
to  put  them  in  a  form  suitable  for  a  reservoir  thermal  simulation  model. 
In  the  experiments,  several  parameters  were  measured,  including  the 
total  volume  flux  entrained  into  the  jet.  For  the  present  purpose,  how¬ 
ever,  it  is  more  convenient  to  fit  an  empirical  equation  to  the  mea¬ 
sured  distribution  of  volume  flux  into  the  jet.  The  parameters  of  the 
equation  are  chosen  to  maintain  consistency  between  the  total  measured 
flux  and  that  computed  from  the  empirical  distribution. 

36.  The  basic  equation  is  (Equation  A56): 


.1/2 


(18) 


where 

q  =  the  entrained  flow  rate  per  unit  depth 
e 

n  =  z/S.  ,  the  nondimensional  height 
£ 

r]  =  parameter  to  account  for  asymmetry  of  entrained  volume 
0  flux  distribution  about  the  jet  axis 


20 


COMPUTE  AND  INCREMENT  ENTRAINED 
VOLUME  FLUX  FOR  EACH  LAYER 


b  =  parameter  to  account  for  width  of  the  entrained 
distribution 

q  =  maximum  of  q  ;  occurs  at  n  =  n 
em  e  o 

The  choice  of  empirical  equations  to  predict  the  parameters  of  Equa¬ 
tion  18  is  discussed  in  Appendix  A  (section.  Entrained  Volume  Flux  Dis¬ 
tribution),  and  the  equations  are  summarized  in  Table  A10.  Comparisons 
of  the  volume  flux  distribution  predicted  by  these  equations  with  the 
actual  measured  distributions  are  shown  in  Figures  All,  A12,  and  A22-A27. 
The  fit  in  all  cases  is  sufficiently  close  for  the  present  modeling  pur¬ 
poses,  and  the  total  entrained  volume  flux  computed  from  these  equations 
is  within  11  percent  of  that  measured.  The  entrained  volume  flux  is 
withdrawn  selectively  from  a  vertical  zone  of  finite  thickness  (see  Fig¬ 
ure  Al).  The  limits  of  this  zone  are  z  and  z  +  W  ,  as  computed 
from  Equations  A43  and  A44  (Table  A10). 

37.  The  entrainment  from  a  zone  whose  bottom  height  is  z^  and 
top  height  is  z^  is: 


=/  qedz 


substituting  Equation  18  into  Equation  19  and  integrating  yields: 

Qe  =  [erf(a2)  -  erfUj)]  (2C 

1/2 

where  q  is  q  /M  ,  CT  is  the  standard  deviation  of  the  assumed 
o  em 

Gaussian  distribution,  Equation  18,  i.e., 
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and  erf  is  the  error  function: 


erf(t)  =  —  f  e  S  ds  (25) 

VF  J 

o 

38.  The  entrainment  distribution  is  computed  in  the  subroutine 
WBJET  in  a  manner  shown  in  the  flowchart  of  Figure  9.  As  the  experimen¬ 
tal  results  apply  to  linear  stratifications,  the  first  step  is  to  lin¬ 
earize  the  ambient  density  profile.  This  is  done  by  a  least  squares  fit 
over  the  depth  corresponding  to  the  top  and  bottom  of  the  jet  outlet, 
and  the  stratification  parameter  £  (see  Appendix  A)  is  computed  using 
this  depth.  The  length  scales  £q  ,  £^  ,  and  £^  are  calculated  from 

Equations  A6 ,  A7 ,  and  A8.  The  vertical  extent  of  the  jet  influence, 
to  ze  +  we  >  calculated  from  Equations  A43  and  A44,  and  the  strati¬ 
fication  parameter  e  is  then  recalculated  over  this  height.  This  pro¬ 
cedure  is  repeated  until  the  calculation  of  £  is  within  1  percent  of 
the  previously  calculated  value.  The  entrained  flow  rate  is  then  com¬ 
puted  for  each  layer  from  Equation  20  where  the  parameters  are  calcu¬ 
lated  from  the  equations  of  Table  A10  and  Equations  21-24. 


Summary  and  Entrainment  Subroutine  Flowcharts 


39.  Flow  parameters  are  first  computed  in  the  subroutine  ENTRAIN. 
The  stratification  parameter  £  is  first  computed  by  linearizing  the 
stratification  profile  over  the  depth  corresponding  to  the  top  and 


bottom  of  the  jet  outlet.  The  length  scales  £ 


Q  ’ 


M 


and  the 


ratios  £q/££  an<*  £^/££  are  then  computed.  As  discussed  in  Appen¬ 
dixes  A  and  B,  jets  for  which  £^/££  <  2  are  strongly  buoyant.  For 
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this  case,  the  subroutine  SBJET  is  called.  If  £„/£  >  2  ,  the  jet  is 

momentum  dominated,  and  the  subroutine  WBJET  (Weakly  Buoyant  JET)  is 
called.  The  flow  diagram  for  ENTRAIN  is  shown  in  Figure  10. 

40.  The  entrainment  distribution  is  computed  in  SBJET  by  a  modi¬ 
fied  version  of  Fan  and  Brooks'  (1969)  integral  model.  The  method  of 
solution  used  in  their  model  is  given  in  Appendix  C;  the  adaptation  to 
WESTEX  is  described  in  paragraphs  25-33.  (The  flow  diagram  for  SBJET  is 
shown  in  Figure  8.) 

41.  The  entrainment  distribution  for  weakly  buoyant  jets  is  com¬ 
puted  in  WBJET  by  use  of  experimental  results  obtained  during  the  pres¬ 
ent  study.  These  results  are  presented  in  Appendix  A,  and  their  adapta¬ 
tion  to  WESTEX  is  described  in  paragraphs  35-38.  (The  flow  diagram  for 
WBJET  is  shown  in  Figure  9.) 

42.  One  final  modification  is  required  to  account  for  the  presence 

of  the  free  surface  or,  as  in  the  case  of  Carters  Lake,  Ga.,  a  nearby 
solid  lower  boundary.  If  +  W^  is  above  the  free  surface,  entrain¬ 
ment  above  is  cut  off;  if  the  location  of  is  below  a  solid  boundary, 

entrainment  below  is  cut  off.  This  assumption  will  generally  reduce  the 
total  flow  entrained  into  the  jet  and  is  used  in  both  SBJET  and  WBJET. 
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Figure  10.  Flowchart  for  subroutine  ENTRAIN 


PART  IV:  APPLICATION  TO  CARTERS  LAKE,  GEORGIA 


43.  The  general-purpose  entrainment  algorithms  presented  herein 
were  incorporated  into  the  reservoir  thermal  simulation  code  (WESTEX). 

The  WESTEX  code  already  contained  the  algorithms  to  simulate  the  effects 
of  PS  operations,  as  discussed  in  Part  II.  To  test  the  capability  of 
the  general  purpose  entrainment  algorithms  within  the  thermal  modeling 
framework,  the  revised  code  was  tested  with  a  site-specific  application. 

A  Corps  of  Engineers  PS  reservoir,  Carters  Lake,  Ga.,  was  selected  for 
the  test  case. 

44.  Carters  Lake  is  a  500-MW  PS  project  located  on  the  Coosawattee 

River  (Figure  11).  At  normal  summer  pool  elevation  (327  m  mean  sea 

6  3 

level  (msl)),  Carters  Lake  contains  473  ><10  m  of  water  with  a  maximum 

depth  of  124  m.  A  reregulation  dam  located  about  2.9  km  below  Carters 

6  3 

Dam  impounds  about  21.6  x  io  m  at  a  maximum  pool  depth  of  approxi¬ 
mately  15.2  m. 

45.  The  maximum  generation  and  pumpback  flow  rates  at  Carters  Dam 

3 

are  611.6  and  213.2  m  /sec,  respectively.  With  these  flows  and  the 
reservoir  depths  and  volumes  noted  above,  the  upper  reservoir  maintains 
a  strong  thermal  stratification,  while  the  lower  reservoir  is  practi¬ 
cally  fully  mixed.  The  heat  budget  of  the  afterbay  (reregulation  pool) 
was  modeled  by  assuming  the  afterbay  was  fully  mixed. 

46.  The  four  penstocks  for  the  four  turbines  are  located  at  invert 
elevation  299  m  at  the  bottom  and  end  of  a  rock-cut  channel  in  the  right 
abutment.  The  entrance  channel  is  approximately  488  m  long  and  a  mini¬ 
mum  of  61  m  wide.  Two  of  the  turbines  are  reversible  pump-turbines  that 
discharge  into  the  upper  reservoir  through  two  penstocks.  The  dimen¬ 
sions  of  each  penstock  intake  are  6.25  m  high  by  4.27  m  wide.  The  length 

scale  relationships  (Equations  A6-A8)  and  ratios  (!L/£  and  £-./£  ) 

y  E  rl  e 

were  computed  for  the  Carters  Lake  pumped  inflows.  During  periods  of 
upper-lake  stratification,  the  PS  inflow  jet  was  found  to  be  within  the 
range  of  applicability  of  the  results  of  Appendix  A  (high  momentum  and 
low  buoyancy,  £^/££  >  2  ;  negligible  effects  of  nozzle  size,  £q/££ 

<  0.5).  During  periods  of  little  or  no  upper-lake  stratification, 


Carters  Lake  project  features 


buoyancy  was  expected  to  be  more  important  (£  /£  <  2)  ,  thus  requiring 

the  use  of  the  integral  buoyant  jet  model,  SBJET. 

47.  The  two  pumpback  outlets  were  modeled  as  one  outlet  with  a 

2  3  2 

discharge  area  of  26.7  m  for  Q  <  113  m  /sec  and  53.4  m  for  Q  >  113 

3  3 

m  /sec  .  The  pumpback  capacity  of  each  pump-turbine  is  about  113  m  /sec. 

The  pumpback  jet  was  expected  to  collapse  before  reaching  the  end  of  the 

rock-cut  channel;  thus,  entrainment  of  ambient  reservoir  water  was  not 

expected  below  the  channel  bottom  elevation.  The  entrainment  algorithms 

were  modified  to  not  allow  any  entrainment  below  the  channel  invert, 

elevation  298  m. 

48.  Hydrological,  meteorological,  and  operational  data  for  1979 
were  used  for  the  simulation.  Vertical  profiles  of  temperature  had  been 
measured  biweekly  at  Carters  Lake  between  May  and  December  1979.  The 
model  simulation  was  started  for  1  January  1979,  with  isothermal  condi¬ 
tions  of  10°  C.  A  time  step  of  a  day  was  used.  Predicted  and  observed 
temperature  profiles  are  compared  in  Figure  12.  In  general,  the  pre¬ 
dicted  profiles  exhibit  too  much  mixing.  It  appeared  that  too  much  en¬ 
trainment  took  place  in  the  simulation. 

49.  Considering  the  constriction  of  the  approach  channel,  it 
seemed  quite  reasonable  that  the  channel  sidewalls  (as  well  as  the  bot¬ 
tom,  as  discussed  above)  could  restrict  the  amount  of  entrainment.  All 
of  the  results  presented  herein  are  based  on  wall  jets  issuing  into  an 
unconfined  fluid. 

50.  The  geometry  shown  in  Figure  11  definitely  does  not  satisfy 
this  criterion,  as  the  jet  is  confined  by  the  sidewalls.  Based  on  the 
experiments  conducted  on  nonbuoyant  jets  in  two-layer  stratification, 
the  bottom  boundary  will  probably  not  affect  the  total  entrainment  much, 
but  the  entrained  flow  distribution  becomes  unsymmetrical ,  with  the  peak 
shifted  toward  the  bottom  (see  Roberts  and  Matthews  1984).  To  compen¬ 
sate  for  this,  the  coefficient  of  0.8  in  Equation  A60  was  gradually 
reduced,  and  the  effect  on  predicted  temperature  profiles  was  investi¬ 
gated.  The  agreement  between  predicted  and  observed  profiles  was  best 
with  a  coefficient  of  zero  rather  than  0.8. 

51.  The  sidewall  effect  of  the  narrow  entrance  channel  is 
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difficult  to  estimate.  It  is  logical  that  the  entrainment  is  reduced. 
The  maximum  entrainment  per  unit  depth  was  reduced  50  percent  by  using  a 
coefficient  of  0.52  rather  than  1.04  in  Equation  A62.  These  two  changes 
to  Equations  A60  and  A62  resulted  in  about  a  60-percent  decrease  in 
total  entrainment.  The  favorable  results  of  these  changes  are  shown  in 
Figure  13. 

52.  The  entrainment  algorithms  functioned  very  satisfactorily  ex¬ 
cept  for  the  site-specific  changes  discussed  above.  Subroutine  SBJET 
was  called  from  simulation  day  1  (1  January  1979)  through  day  58  (27  Feb 
ruary) .  After  day  58,  subroutine  WBJET  was  called  every  day  through  the 
rest  of  the  simulation  (end  of  simulation  was  day  341),  except  for  a 
couple  of  days  in  the  spring.  This  indicates  that  at  Carters  Lake  for 
1979,  the  stratification  was  strong  enough  by  early  March  to  cause  the 
jet  to  be  weakly  buoyant. 
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Figure  13.  Carters  Lake  predicted  and  observed 
temperature  profiles  with  site-specific  modifi 
cations  to  entrainment  equations 


PART  V:  SUMMARY  AND  RECOMMENDATIONS 


53.  Pumped-storage  inflows  can  have  a  significant  impact  on  reser¬ 
voir  thermal  stratification.  In  the  case  of  Carters  Lake,  the  epilim- 
nion  and  metalimnion  were  deepened  to  a  depth  of  about  30  m,  compared 
with  about  10  m  which  is  typical  of  most  lakes  in  the  region.  This  mix¬ 
ing  and  deepening  of  the  thermocline  have  been  accompanied  by  a  substan¬ 
tial  increase  in  the  total  heat  in  the  reservoir. 

54.  The  entrainment  process  must  be  modeled  in  order  to  simulate 
the  mixing  that  occurs  with  PS  inflows.  To  illustrate  this,  a  simula¬ 
tion  was  conducted  for  Carters  Lake  without  any  entrainment.  The  re¬ 
sults  of  this  simulation  are  shown  in  Figure  14,  along  with  the  observed 
data  and  the  predictions  with  entrainment.  The  differences  in  the  pre¬ 
dictions  without  entrainment  and  those  with  entrainment  (and  the  ob¬ 
served)  indicate  that  significant  mixing  occurs  due  to  the  entrainment 
process  and  that  this  process  must  be  modeled  to  provide  realistic 
predictions . 

55.  The  research  described  herein  provided  a  general-purpose  en¬ 
trainment  algorithm.  The  algorithm  should  be  applicable  to  most  PS 
inflow  jets,  provided  that  the  inflow  is  a  horizontal  or  inclined 

(0  <  45  deg)  wall  jet  issuing  into  an  unconfined  fluid.  Lateral  and 
vertical  restrictions  in  the  fluid  surrounding  the  jet  impose  site- 
specific  considerations  as  with  the  Carters  Lake  application  presented 
here.  The  Carters  Lake  intake  channel  geometry  was  restrictive,  but  it 
did  not  preclude  the  use  of  the  algorithm,  although  it  did  require  three 
modifications.  One  modification  was  made  to  prevent  any  entrainment  be¬ 
low  the  channel  bottom  (see  Part  IV).  The  other  two  were  made  to  reduce 
the  total  amount  of  entrainment  and  to  modify  the  location  of  peak  en¬ 
trained  flow,  as  these  parameters  were  affected  by  sidewall  and  bottom 
boundaries,  respectively.  Permanent  coding  changes  were  implemented  to 
account  for  the  bottom  effects.  The  effects  of  sidewall  restrictions 
must  be  handled  on  a  site-specific  basis. 

56.  If  sidewall  interference  is  questionable,  then  the  collapse 
length  can  be  used  as  a  first  approximation  for  the  nonrestrictive  zone. 
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The  collapse  length  can  be  estimated  by  3j0£  (see  Table  Bl).  If 

either  sidewall  is  closer  than  about  x^/3  »  or  ,  from  the  jet 

center  line,  lateral  restriction  will  probably  occur.  If  lateral  re¬ 
striction  is  anticipated,  then  the  coefficient  in  Equation  A62  must  be 
reduced.  The  amount  of  reduction  must  be  calibrated  with  observed  pro¬ 
files.  If  observed  profiles  do  not  exist,  it  would  be  advisable  to 
examine  the  sensitivity  of  this  coefficient  and  to  try  to  place  some 
bounds  on  the  predictions.  An  undistorted,  near-field  physical  model 
could  be  used  to  measure  the  entrainment  as  influenced  by  local  geometry 
(e.g.,  sidewalls).  The  physical  model  results  would  be  used  to  adjust 
coefficients  in  the  entrainment  algorithm. 
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APPENDIX  A:  RESULTS  OF  EXPERIMENTS  ON  WEAKLY  BUOYANT  JETS 


Introduction 


1.  Investigation  of  the  literature  on  weakly  or  neutrally  buoyant 
jets  discharging  into  a  stratified  ambient  revealed  little  information. 
Experiments  were  therefore  conducted  to  study  this  situation.  As  even 
the  restricted  idealized  problem,  shown  in  Figure  Al,  represents  a  large 
number  of  experiments  when  all  possible  combinations  of  jet  geometries, 
angles,  density  stratifications,  and  buoyancies  are  considered,  emphasis 
was  placed  on  the  following  categories: 

a.  Series  1.  Horizontal  neutral  and  slightly  buoyant  dis¬ 
charges  into  a  linear  stratification.  The  nozzle  size 
was  also  varied  to  investigate  the  effects  of  this  param¬ 
eter.  Results  are  presented  in  paragraphs  20-22. 

b.  Series  2.  Inclined  jets  discharging  at  22.5  or  45  deg 
to  the  horizontal  into  a  linear  stratification.  Buoy¬ 
ancies  were  neutral,  slightly  positive,  or  slightly  nega¬ 
tive.  The  nozzle  size  was  maintained  small  to  eliminate 
effects  of  this  parameter.  Results  are  presented  in 
paragraphs  23-26. 

c.  Series  3.  Horizontal  neutrally  buoyant  jets  discharging 
into  the  lower  layer  of  a  two-layer  fluid.  The  upper 
layer  is  well  mixed  and  separated  from  the  linearly  stra¬ 
tified  lower  layer  by  a  finite  density  jump.  The  depth 
of  discharge  below  the  interface  was  varied,  but  the 
nozzle  size  was  kept  small.  Results  are  presented  in 
paragraphs  27-29. 

2.  The  dynamics  of  nonbuoyant  jets  in  uniform  density  fluids  have 
been  extensively  studied  (i.e.,  the  review  in  Fischer  et  al.  1979)*  and 
their  gross  characteristics  are  well  known.  Such  a  jet  grows  linearly 
with  distance  from  the  source,  and  its  total  volume  flux  continuously 
grows  with  distance.  Schneider  (1980)  performed  experiments  on  hori¬ 
zontal  nonbuoyant  jets  discharging  into  the  lower  layer  of  a  stagnant 
fluid  consisting  of  two  homogeneous  layers  of  different  densities.  He 
measured  the  entrainment  rate  across  the  interface  and  showed  that  this 

*  Bibliographic  information  is  given  in  the  References  section  of  the 
main  text,  page  36. 


rate  is  a  function  of  a  form  of  the  Richardson  number  based  on  the  depth 
of  the  discharge  below  the  thermocline,  the  density  jump  across  the 
thermocline,  and  the  jet  momentum  flux.  Jirka  et  al.  (1977)  reported 
experiments  on  jet  discharges  from  ocean  thermal  energy  conversion 
plants  into  the  upper  layer,  which,  although  similar  to  the  present 
situation,  cannot  be  compared  directly  due  to  their  emphasis  on  flows  in 
the  uniform  upper  layer. 

3.  When  a  jet  is  discharged  into  a  continuously  stratified  medium, 
its  behavior  changes  completely  from  one  discharged  into  a  homogeneous 
fluid.  Near  to  the  source,  the  jet  growth  is  unaffected  by  stratifica¬ 
tion,  but  further  away  the  initial  jet  velocity  decays  and  the  width 
grows  until  the  local  gradient  Richardson  number  becomes  large  enough 
that  the  stratification  removes  the  turbulent  kinetic  energy  and  the  jet 
collapses  vertically,  causing  entrainment  to  essentially  cease.  The  en¬ 
trained  flow  is  selective  due  to  the  ambient  stratification.  Following 
collapse,  the  inflow  moves  upstream  as  a  tonguelike  density  current. 

4.  In  the  following  section  a  general  dimensional  analysis  of  the 
problem  is  presented  in  terms  of  length  scales.  The  experimental  pro¬ 
cedure  is  discussed  in  paragraphs  14-19,  and  the  results  are  given  in 
paragraphs  20-29.  A  brief  discussion  of  the  results  and  equations  to 
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describe  them  are  presented  in  paragraphs  30-37.  The  distribution  of 
entrained  volume  flux  for  all  the  experiments  and  equations  suitable 
for  inclusion  into  the  pumped-storage  (PS)  model  are  given  in  paragraphs 
38-43.  Concluding  remarks  are  given  in  paragraphs  44-45. 

Analysis 

5.  The  problem  under  consideration  is  shown  in  Figure  Al.  It  is 
most  easily  analyzed  and  discussed  by  means  of  the  characteristic  length 
scales  involved,  following  Fischer  et  al.  (1979).  Assuming  the  jet  to 
be  fully  turbulent  so  that  Reynolds  number  effects  are  negligible,  the 
jet  can  be  characterized  by  its  kinematic  fluxes  of  volume  Q  ,  momentum 
M  ,  and  buoyancy  B  ,  where: 

Q  =  uA  (Al) 


M  =  u2A 


B  =  uA 


where 


(^) 


u  =  jet  nozzle  velocity 
A  =  nozzle  area 

=  ambient  density  at  the  discharge  depth 
Pq  =  jet  density 
g  =  acceleration  due  to  gravity 

6.  By  making  use  of  the  Boussinesq  assumption  (Turner  1973),  one 
can  characterize  the  ambient  stratification  by  g'  and  £  where 
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and  Ap  is  the  density  jump  across  the  interface.  The  square  root  of 
£  is  the  Brunt-Vaisala,  or  buoyancy,  frequency. 

7.  Many  characteristic  length  scales  can  be  formed  from  the  six 
parameters  involved.  This  study  is  concerned  with  situations  in  which 
the  momentum  flux  is  always  important  and  the  volume  and  buoyancy  fluxes 
less  so.  Five  independent  length  scales  can  be  formed,  the  most  con¬ 
venient  of  which  are: 
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8.  The  problem  is  characterized  by  the  relative  magnitudes  of  h 
(distance  from  elevation  of  jet  source  to  density  jump)  and  these  five 


length  scales,  each  of  which  has  a  physical  significance.  The  term  SL 


is  the  square  root  of  the  nozzle  area;  it  characterizes  the  distance 
from  the  origin  for  which  the  nozzle  geometry  affects  the  flow  field. 


For  distances  from  the  origin  much  greater  than  £q  ,  the  properties  of 


the  jet  depend  only  on  M  and  distance  and  are  independent  of  the 


nozzle  geometry  (Fischer  et  al.  1979,  p  324).  The  term  JL^  is  the 


characteristic  distance  from  the  nozzle  over  which  the  source  momentum 
flux  is  important  relative  to  the  source  buoyancy  flux.  In  an  unstra¬ 


tified  environment,  for  distances  from  the  origin  much  less  than  SL 


M 


the  momentum  flux  dominates  and  the  flow  is  jetlike;  for  distances  much 


greater  than  ,  the  buoyancy  flux  dominates  and  the  flow  is  plume¬ 


like.  For  distances  of  the  order  of 


£„  ,  both  momentum  and  buoyancy 
M 


fluxes  are  important.  The  term  ££  characterizes  the  strength  of  the 


jet  relative  to  the  ambient  linear  stratification.  The  same  length  scale 


arises  in  vertically  discharged  jets  in  a  linearly  stratified  ambient 
(Fischer  et  al.  1979),  in  which  case  the  terminal  rise  height  is  directly 
proportional  to  £  The  term  £  characterizes  the  relative  strengths 

£  U 

of  the  jet  momentum  to  the  interfacial  density  jump. 

9.  There  are  several  parameters  of  interest  for  the  present  study 

(see  Figure  Al).  These  include  the  collapse  length  X  of  the  jet;  the 

m 

location  of  the  neutrally  buoyant  collapsed  layer  Z.  ;  the  collapsed 

L 

layer  thickness  Wr  ;  the  location  Z  and  width  W  of  the  entrained 
L  e  e 

flow;  the  mixed  layer  thickness  W  ,  defined  as  the  thickness  of  the 
ambient  fluid  where  the  density  profile  is  altered  by  the  jet  after  the 
inflow  ceases;  and  Z  ,  the  location  of  the  bottom  of  this  layer.  It 
is  also  desirable  to  know  whether  fluid  is  mixed  into  the  jet  from  the 
upper  layer  in  the  two-layer  case.  This  will  be  termed  stability,  where 
stable  signifies  no  entrainment  from  the  upper  layer  and  unstable  sig¬ 
nifies  that  entrainment  from  the  upper  layer  does  occur. 

10.  The  dependent  variables  can  be  expressed  as: 

V  WL’  V  ze’  Wm ’  Zm’  ZL  =  £(Q’  M’  B’  h’  £>  8' ’0)  (A10) 

and 

stability  =  f (Q,  M,  B,  h,  e,  g’ ,  6)  (All) 


Following  dimensional  analysis,  Equations  A10  and  All  can  be  written  in 
terms  of  the  length  scales  as: 


and 


(A12) 


stability 


(A13) 


11.  The  total  volume  flux  in  the  jet  becomes  finite  following  col¬ 
lapse  and  does  not  increase  continually  with  distance  as  in  the  unstra¬ 
tified  case.  Defining  this  total  volume  flux  as  p  ,  we  have: 


A5 


M  =  f (Q,  M,  B,  h,  e,  g\  0) 


which,  following  dimensional  analysis  becomes: 


J£1/4  =  f  (!2  !l  0 

b3/4  -  W  lc-  V  V 


(A14) 


(A15) 


The  term  is  chosen  as  the  normalizing  length  scale  as  this  will 

always  be  important  for  the  present  class  of  problems. 

12.  The  entrainment  distribution  can  be  derived  by  similar  argu¬ 
ments.  Suppose  that  is  the  entrained  volume  flow  rate  per  unit 

depth,  so  that  the  total  entrained  flow  rate  is: 


Z  +W 
e  f  e 


■I 


q  dz 
e 


(A16) 


then: 


qe  =  f(Z,  M,  B,  Q,  h,  e ,  g* ,  0) 


(A17) 


which,  following  dimensional  analysis  and  expression  in  terms  of  length 


scales,  becomes: 


_^  =  ffz_!2!Mh_h_ 

mi/2  i  r  r  i’  r 


and  the  maximum  value  q  is: 

em 


*em  _  f  (!Q  h_  h_ 

1/2'  \  *’  V  £’  £’ 

M  \  e  e  e  z  t 


(A18) 


(A19) 


13.  Equations  A12,  A13,  A15,  A18,  and  A19  are  based  on  dimensional 
analysis  only  and  could  be  expressed  in  terms  of  the  fundamental  vari¬ 
ables  Q  ,  M  ,  B  ,  h  ,  e  ,  and  g'  .  It  is  more  convenient  to  ex¬ 
press  them  in  terms  of  the  length  scales,  however,  due  to  their  physical 
significance  and  aid  to  interpretation. 


A6 


Experimental  Procedure 

14.  A  sketch  of  the  experimental  configuration  is  shown  in  Fig¬ 
ure  A2.  The  experiments  were  performed  in  a  Plexiglas-walled  tank  8  ft 
(2.4  m)  long,  4  ft  (1.2  m)  wide,  and  2  ft  (0.6  m)  deep.  The  discharge 
was  from  a  round  nozzle.  The  stratification  was  obtained  with  salt  by 
means  of  a  two-tank  filling  method  similar  to  that  used  by  Wright 
(1977a).  The  fluid  was  introduced  to  the  tank  by  a  floating  spreader. 
Salt  concentrations,  and  hence  density  profiles,  were  measured  by  a  two- 
wire  conductivity  probe  mounted  on  a  vertically  traversing  probe  on  an 
instrument  carriage.  A  Sanborn  recorder  was  used  to  monitor  the  conduc¬ 
tivity  probe;  the  probe  and  circuitry  were  very  similar  to  those  used  by 
Wright  (1977a).  The  jet  fluid  was  dyed  to  make  it  visible,  and  photo¬ 
graphs  of  the  flow  field  were  obtained  by  an  overhead  motor-driven  35mm 
camera.  A  mirror  placed  at  the  side  of  the  tank  enabled  simultaneous 
side  and  overhead  views  to  be  obtained  in  the  same  photograph. 

15.  Due  to  the  small  velocities  involved,  it  is  quite  difficult  to 
measure  the  entrained  volume  flux.  It  was  done  by  taking  advantage  of 
the  entrainment-induced  flow  field  in  the  tank.  Prior  to  commencement 

of  the  experiment,  potassium  permanganate  dye  crystals  were  dropped  into 
the  tank  along  a  straight  line  perpendicular  to  the  jet  axis,  and  beyond 
the  point  at  which  the  jet  would  collapse.  The  subsequent  movement  of 
the  dye  streaks  was  photographed,  and  the  entrained  flow  was  computed. 

16.  As  not  all  of  the  entrainment-induced  flow  is  actually  en¬ 
trained  into  the  jet,  it  was  necessary  to  measure  independently  the 

width  of  the  entrained  flow  W  and  its  location  Z  .  This  was  done 

e  e 

by  placing  the  conductivity  probe  within  the  jet,  first,  near  its  top 
and,  second,  near  its  bottom  to  record  the  salinity  fluctuations.  As 
material  entrained  into  the  jet  is  not  immediately  mixed  (Roshko  1976), 
the  maxima  and  minima  of  the  salinity  fluctuations  at  any  level  corre¬ 
spond  to  the  upper  and  lower  depths  from  which  the  fluid  was  entrained. 

17.  Salinity  extremes  were  calculated  and,  from  the  measured  den¬ 
sity  profile,  the  heights  corresponding  to  these  concentrations  were  com¬ 


puted.  This  gives  the  limits  of  the  entrainment  width  and  its  location, 


entrained  flow  rate  per  unit  horizontal  width,  and  the  total  entrained 

flow  Q  is 
e 

Q  =  slope  x  w  (A21) 

6 

where  w  is  the  horizontal  width  of  the  entrained  flow  field  that  was 
estimated  from  the  overhead  photographs.  The  flow  field  was  found  to  be 
very  uniform  across  the  tank  width  beyond  the  jet  boundaries.  Finally, 
the  total  volume  flux  p  in  the  jet  is  computed  by  conservation  of 
volume  as: 

p  =  Qfi  +  Q  (A22) 

19.  The  thickness  of  the  collapsed  layer  and  the  approximate 

length  were  measured  from  the  photographs. 

Results 

Horizontal  buoyant  and  nonbuoyant 
jets  into  linear  stratification 

20.  For  these  experiments  there  is  no  interfacial  density  jump; 
thus,  SL ^  (Equation  A9)  and  h  are  not  parameters.  Experiments  were 
performed  for  the  following  parameter  ranges: 


0.08  <  ^  <  1 

e 

and 

£m 

r  =  2  or  oo 
£ 

For  these  experiments  W  and  Z  ,  the  mixed  layer  thickness  and  loca 

ID  ID 

tion,  were  not  measured.  For  0=0  deg  ,  Equations  A12,  A15,  A18,  and 
A19  become: 


X  WT  zr  w  z 

m  L  L  e  e 
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(A23) 


(A2A) 


(A25) 


(A26) 


21.  Side-view  photographs  of  neutral  (2/2  =  “)  and  slightly 

rl  £ 

buoyant  (2^/2£  =  2)  jets  are  shown  in  Figures  A3  and  AA,  and  a  simul¬ 
taneous  overhead  and  side  view  is  shown  in  Figure  A5  for  a  neutrally 
buoyant  case.  Many  of  the  characteristics  of  all  the  jets  studied  are 
apparent  in  these  photographs,  including  the  initial  unimpeded  jet 
growth,  the  sudden  vertical  collapse  and  sideways  spreading,  the  tongue¬ 
like  horizontal  intrusion,  and  the  entrainment  dye  streak.  Experimental 
parameters  for  this  series  of  experiments  are  given  in  Tables  A1  and  A2. 

22.  Results  for  collapsed  layer  thickness  ,  entrainment  width 

,  entrainment  location  Zg  ,  and  rise  height  Z^  for  the  buoyant 

case  are  shown  in  Figures  A6-A9,  expressed  in  the  form  of  Equation  A23. 
The  total  entrained  volume  flux,  Equation  A2A,  is  shown  in  Figure  A10. 
The  distribution  of  entrained  volume  flux  is  shown  plotted  in  the  form 
of  Equation  A25  in  Figure  All  for  nonbuoyant  jets  =  00  )  and  fig¬ 
ure  A12  for  buoyant  jets  (2„/2  =  2).  The  implications  of  these  results 

M  £ 

are  discussed  in  paragraphs  30-37. 

Inclined  buoyant  and  nonbuoyant 
jets  into  linear  stratification 

23.  The  results  of  the  previous  section  show  the  effects  of  the 
nozzle  size  2^  to  be  negligible  for  2^/2^  <  ~0-2  •  Even  for  larger 
values  of  2/2  ,  the  results  for  W  ,  W  ,  Z  ,  and  Z  (Figures 

y  £  ii  6  C  la 

A6  through  A9)  are  essentially  independent  of  this  parameter.  Only  the 
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Figure  A8.  Experimental  results,  bottom  location  of  entrainment  layer 
caused  by  horizontal  neutrally  and  slightly  buoyant  jets  in  linear 

stratification 
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Figure  A 10.  Experimental  results  for  entrained  volume  flux  of  hori 
zontal  neutral  and  slightly  buoyant  jets  in  linear  stratification 


-1.8 


Figure  A12.  Distribution  of  entrained  volume  flux  for  horizontal 

slightly  buoyant  OL./JZ  =  2)  jets  in  linear  stratification 
M  £ 
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results  for  volume  flux  (Figure  A10)  show  a  dependence  on  £^  ,  and  these 


also  approach  a  constant  value  for  small 


This  result  is  ex¬ 


pected,  since  for  distances  from  the  origin  much  greater  than  £q  ,  the 
properties  of  a  jet  are  dominated  by  the  momentum  flux  M  and  become 
independent  of  the  nozzle  size  (Fischer  et  al.  1979,  Chap.  9). 

24.  Therefore,  this  series  of  experiments  was  conducted  for  a 

small  nozzle  size  only,  i.e.,  £ /£  <  0.1  .  Only  the  0.25-in.  (6.3-mm) 

Q  £ 

nozzle  was  used.  The  angles  of  inclination  0  were  22.5  deg  and  45  deg, 

with  buoyancies  corresponding  to  £„/£  =  00  (nonbuoyant)  and  2  (posi- 

n  £ 

tively  and  negatively  buoyant).  The  term  Jd  / £  =  -2  denotes  a  nega- 

tively  buoyant  jet.  The  experimental  parameters  for  all  experiments  are 
shown  in  Tables  A3  and  A4. 

25.  Photographs  of  typical  jets  are  shown  in  Figures  A13  and  A14. 
Some  of  the  jet  characteristics,  in  particular  the  jet  collapse  and 

the  rapid  spreading  after  collapse,  are  similar  to  those  for  horizontal 
jets,  as  discussed  in  the  previous  section. 

26.  For  this  class  of  jets  £ is  again  not  a  parameter,  as  there 
is  no  interfacial  density  jump.  As  previously  discussed,  the  nozzle 
size  is  not  a  parameter,  and  so  Equations  A12,  A15,  and  A18  become: 


and 


(A29) 


The  experimental  results  for  these  parameters  are  plotted  in  the  form  of 
Equations  A27,  A28,  and  A29  in  Figures  A15  through  A27.  These  results 
are  discussed  further  in  paragraphs  30-37. 


Figure  A13.  Side-view  photograph  of  inclined  jet  (0  =  45  deg) 
into  linear  stratification 


Figure  A14.  Simultaneous  overhead  and  side  view  (in  mirror  at  top) 
of  inclined  jet  (6  =  45  deg)  into  linear  stratification 
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Figure  A17.  Experimental  results  for  thickness  of 
entrainment  layer  of  inclined  jets  of  slightly 
negative  (£M/*£  =  -2),  neutral  =  ®)  >  and 

slightly  positive  buoyancy  =  2)  into 

linear  stratification 


<■**]  tNJ 


-2.00  A 

oo 

+  2  00 


o  a 


Figure  A21.  Experimental  results  for  entrained  volume 
flux  of  inclined  jets  of  slightly  negative  (£mM£  =  -2)» 

neutral  (£.,/.£  =  “) ,  and  slightly  positive  buoyancy 

n  £ 

(£. /£  =  2)  into  linear  stratification 

M  £ 


re  A25.  Distribution  of  entrained  volume  flux  into  inclined  jet, 
0  =  45  deg,  slightly  buoyant  =  2) 


A30 


re  A26.  Distribution  of  entrained  volume  flux  into  inclined  jet 
0  =  45  deg,  neutrally  buoyant  (£^/££  =  00 ) 
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Horizontal  nonbuoyant  jets 
into  two-layer  stratification 

27.  The  form  of  the  stratification  for  these  experiments  is  as 
shown  in  Figure  Al.  A  well-mixed  upper  layer  overlies  a  linearly 
stratified  lower  layer.  The  layers  are  separated  by  a  density  jump  of 
magnitude  Ap  .  Only  horizontal  nonbuoyant  jets  were  considered.  As  in 
the  previous  section,  the  0.25-in.  (6.3-mm)  nozzle  was  used  and  £  /£ 

<  ~0.1  ,  so  that  effects  of  £^  could  be  neglected.  Equations  A12, 

A13,  and  A15  then  become: 


X  WT  ZT  W  Z  W  2  /  .  \ 

m  L  L  e  e  m  tn  _  -In  h  \ 

£’£’£’  £’£’£’£  “  M  £  ’  £A  J 
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stability 


(A30) 


(A31) 


je1/4  h_  h_ 

H3/4  V  *A 


(A32) 


Equation  A18  becomes: 


=  /z_  h_  h_\ 

mi/2  f\vVV 


(A33) 


28.  Parameters  for  all  experiments  are  shown  in  Tables  A5  and  A6. 
Photographs  of  typical  experiments  are  shown  in  Figures  A28  and  A29.  In 
Figure  A28,  the  jet  is  submerged  well  below  the  density  jump,  and  the 
appearance  of  the  jet  is  very  similar  to  that  in  a  continuous  linear 
stratification,  Figure  A3.  The  jet  in  Figure  A29  is  closer  t  >  the  inter¬ 
face,  the  effect  of  which  is  to  sharply  suppress  the  turbulence  there, 
causing  a  flattened  top  to  the  jet. 
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Figure  A28.  Simultaneous  overhead  and  side  view  (in  mirror 
at  top)  of  horizontal  nonbuoyant  jet  in  two-layer  stratifi 
cation.  Jet  is  well  below  interface. 


Figure  A29.  Simultaneous  overhead  and  side  view  (in  mirror 
at  top)  of  horizontal  nonbuoyant  jet  in  two-layer  stratifi 
cation.  Jet  is  near  interface. 


29.  Experimental  results,  presented  in  the  form  of  Equations  A30, 
A31,  and  A32,  are  given  in  Figures  A30  through  A36.  The  entrained  flow 
distribution  is  shown  in  Figure  A37  for  jets  far  from  the  interface 

(h/£  >  0.7)  and  in  Figure  A38  for  jets  near  the  interface  (h/£  <  0.7). 

£  £ 

Discussion 


30.  A  great  deal  of  data  were  collected  during  the  course  of  the 
three  sets  of  experiments.  In  this  section,  only  the  results  pertinent 
to  the  modeling  of  PS  reservoirs  will  be  discussed.  Further  discussion 
can  be  found  in  a  number  of  journal  articles  prepared  on  these  results 
(Roberts  and  Matthews  1984,  for  example). 

31.  Discussion  of  the  analysis  of  the  effects  of  inflow  angle, 
source  volume,  momentum,  and  buoyancy  fluxes,  and  ambient  stratifica¬ 
tion  parameters  by  dimensional  and  length  scale  arguments  was  presented 
in  paragraphs  5-13.  Perhaps  the  most  important  result  is  that  the  ex¬ 
periments  confirmed  these  results.  The  effect  of  each  parameter  is 
summarized  below,  and  then  limited  equations  are  presented  to  summarize 
the  results.  The  discussion  applies  only  to  the  range  for  which  experi¬ 
ments  were  conducted;  that  is,  high  momentum,  low  buoyancy,  for  which 

SL  /SL  >  2  :  8.  /SL  <  ~1  . 

M  £  -  ’  Q '  z 

32.  The  nozzle  size  does  not  play  a  significant  role,  as  ex¬ 
pected.  For  volume  flux,  the  results  are  independent  of  for  £n/£ 

<0.2  (Figure  A10).  The  effect  on  collapse  width,  entrainment  width, 
entrainment  layer  location,  buoyant  jet  rise  height,  and  entrained  volume 
flux  distribution  is  even  less,  the  results  being  independent  of  £q 

for  £  /£  <  ~0.5  ,  as  shown  in  Figures  A6  through  A9.  For  this  reason, 

Q  £ 

the  effect  of  £^  was  neglected  in  the  subsequent  experiments,  series  2 

and  3,  where  £  /£  was  kept  less  than  0.1. 

Q  £ 

33.  The  dominant  parameters  for  all  experiments  were  the  momentum 
flux  M  and  the  stratification  parameter  £  .  The  collapsed  layer 
thickness  (Figures  A6,  A15,  and  A30)  depended  only  on  M  and  £  .  Other 
parameters  are  also  primarily  controlled  by  M  and  £  ,  although  their 
values  are  slightly  modified  by  the  inclination  angle  and  buoyancy.  For 
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Figure  A3I .  Experimental  results  for  entrained  layer  thickness  of 
horizontal  nonbuoyant  jet  into  two-layer  stratification 


A36 


Figure  A36.  Stability  of  horizontal  nonbuoyant  jet  into 
two-layer  stratification.  Stable  signifies  no  entrainment 
through  the  density  interface. 


example,  the  entrainment  width  and  total  volume  flux  are  independent  of 
buoyancy,  but  vary  with  the  angle  of  the  PS  inflow  (Figures  A7,  A10, 

A17 ,  A21 ,  A31 ,  A35) . 

34.  The  effect  of  the  buoyancy  flux  B  is  primarily  on  the  rise 
height  and  the  location  of  the  entrainment  layer  bottom  (Figures  A8, 

A16,  and  A18) .  The  effect  of  B  on  all  other  parameters  can  be  neg¬ 
lected  for  -  2  •  The  rise  height  and  entrainment  layer  bottom 

are  also  affected  by  the  angle  of  inflow  6  ,  as  shown  in  Figures  Al6 
and  A18.  The  thickness  of  the  entrainment  layer  Wg  is  also  affected 
by  0  ,  becoming  thicker  as  0  increases  (Figure  A17).  The  total  vol¬ 
ume  flux  also  increases  with  0  (Figure  A21),  although  the  dependence 
is  weak. 

35.  The  presence  of  the  interfacial  density  jump  in  the  two-layer 
experiments  (Figures  A30-A35)  does  not  affect  the  gross  flow  parameters 
significantly.  The  collapsed  layer  thickness  is  apparently  reduced  as 
the  jet  is  moved  near  to  the  interface  (Figure  A30).  One  interesting 
effect  is  that  the  maximum  value  of  the  entrained  volume  flux  distribu¬ 
tion  (Figures  A37  and  A38)  does  appear  to  be  substantially  increased  over 
that  in  a  continuous  linear  stratification.  The  reason  for  this  is  not 
presently  known,  but  it  may  be  caused  by  differing  internal  circulations 
and  wave  patterns.  The  scatter  of  the  total  volume  flux  results  (Fig¬ 
ure  A35)  makes  it  difficult  to  deduce  definite  conclusions  on  the  effects 
of  the  interface  on  the  entrained  flow  distribution. 

36.  With  these  observations,  it  is  now  possible  to  deduce  limiting 
formulas  for  the  equations  presented  in  paragraphs  20-29.  For  example, 
Equation  A12  for  the  entrainment  layer  thickness  of  a  horizontal  jet  in 

a  linear  stratification  is: 


For  £q/££  <<  1  ,  and  neutrally  buoyant,  =  00  *  Equation  A34 

becomes : 


f  =  C1  <A35) 

e 

where  is  an  experimental  constant.  The  results  of  Figure  A7  sug¬ 

gest  a  value  of  1.40  for  .  Figure  A7  also  shows  the  value  of 
to  apply  to  slightly  buoyant  jets,  £„/£  =  2  ,  and  to  be  valid  for 

values  of  £  /£  up  to  about  0.5.  Therefore,  Equation  A35  can  be 
written: 


W„  £m  £n 

f  =  1.40  for  ^  >  2  ,  ^  <  0.5 
£  £  £ 


(A37) 


37.  Similar  reasoning  (see  paragraph  36)  can  be  applied  to  the 
other  results  given.  The  resulting  formulas  (Equations  36-50)  are  pre¬ 
sented  in  Tables  A7-A9,  along  with  their  range  of  limitation,  and  are 
plotted  on  the  graphs  where  appropriate. 


Entrained  Volume  Flux  Distribution 


38.  The  distribution  of  entrained  volume  flux  with  height  is  given 
in  dimensionless  form  by  Equations  A18  and  A19,  and  plotted  in  the  pre¬ 
vious  sections.  For  the  purpose  of  using  these  results  for  numerical 
modeling  of  entrainment  in  PS  reservoirs,  it  is  convenient  to  use  empir¬ 
ical  curves  to  describe  this  entrainment  distribution.  Consider  first 
the  horizontal,  neutrally  buoyant  jet  in  linear  stratification.  Equa¬ 
tions  A18  and  A19  for  this  case  become: 


and 


(A51 ) 


(A52) 


39.  Following 
can  be  neglected,  if 


the  arguments  of  paragraphs  30-37,  the  effect 
small.  Thus,  Equations  A51  and  A52  become: 


Y*e 


and 


(A53) 


(A54) 


where  C is  an  experimental  constant.  The  distribution  expressed  by 
Equation  A53  is  shown  in  Figure  All,  which  suggests  a  value  of  1.04  for 
•  We  will  now  assume  that  the  distribution  of  Equation  A53  can  be 
expressed  approximately  as  Gaussian,  i.e.: 


em 


.1/2 


(A55) 


1/2 

where  q^/M  Is  given  by  Equation  A54,  q  =  Z/££  ,  and  b  (see 
Equation  A61)  is  characteristic  of  the  entrainment  width.  Equation  A55 
is  shown  plotted  in  Figure  All  between  the  limits  Zg  to  Zg  +  Wg  with 
a  value  of  b  of  0.65,  and  gives  a  good  fit  to  the  data. 

40.  To  extend  this  distribution  to  the  other  cases,  we  also  assume 
a  Gaussian  distribution  not  centered  on  the  origin,  i.e.: 


q 

e  _  em 

1/2  =  njl/2 


(A56) 


The  task  is  to  find  the  values  of 
cases . 


q  /M 

em 


1/2 


n  ,  and  b  for  all 
o 


41.  Based  on  the  previous  results,  the  rise  height,  represented  by 

1/2 

;  the  entrainment  width,  represented  by  b  ;  and  q=_/M 


size 


;  and 
and  6 

being  of  secondary  importance.  That  is: 


expected  to  be  primarily  dependent  on 


would  be 

'em' 

with  the  nozzle 


A44 


42.  Equation  A56  must  also  satisfy  the  earlier  constraints  on  the 
total  volume  flux,  jj  ,  when  integrated  between  the  lower  and  upper 
limits  of  the  entrained  volume  flux  distribution.  That  is: 


Z  +W 
e  e 

Qe  =  J  qe  dz 

Z 

e 


(A58) 


On  substituting  Equations  A8,  A56,  A6,  and  A22  into  Equation  A58,  we 
obtain,  after  considerable  manipulation: 


«3/4  2H1/2 


z/2  +  w  /«  -  n 


b  er 


-  erf 


Z  /£ 
e  £ 


*  / 


(A59) 


which  must  be  consistent  with  the  earlier  obtained  values  of  total 

volume  flux.  The  empirical  formulas  fitted  to  Equation  A57,  along  with 

those  used  for  W  and  Z  are  shown  in  Table  A10. 
e  e 

43.  The  assumed  entrained  flow  distribution  Equation  A56,  with 
parameters  computed  from  the  equations  given  in  Table  A10,  is  plotted 
along  with  the  experimental  results  in  Figures  A11-A12  and  A22-A27.  The 
assumed  distribution  gives  a  good  fit  in  all  cases.  As  a  further  check, 
the  total  volume  flux  was  computed  from  Equation  A59  and  compared  to  the 
experimental  results  of  Figures  A10  and  A21.  The  values  computed  in 
this  way  were  from  7  to  11  percent  higher  than  those  measured  directly. 
Considering  the  experimental  errors  involved,  this  is  an  acceptable  fit. 


•  •-V'l 


Concluding  Remarks 

44.  A  great  deal  of  experimental  data  were  collected  during  this 

study,  as  summarized  in  this  appendix.  The  results  are  presented  in  a 

form  suitable  for  direct  inclusion  into  a  PS  reservoir  model.  Of  these 

results,  many  are  of  a  fundamental  nature,  and  the  constants  obtained 

are  fundamental.  These  include  the  constants  for  nonbuoyant  horizontal 

jets  into  linear  stratification,  with  SL  / SL  «  1  ,  for  example:  Equa- 

Q  £ 

tions  A36,  A37,  A38,  A41,  and  A42:  and  for  buoyant  jets  with  &*./&  =  2 

n  £ 

Equations  A39  and  A40. 

45.  The  extension  of  the  nonbuoyant  results  to  buoyant  jets  and 
the  other  empirical  formulas  for  inclined  jets  shown  in  Tables  A7-A9  are 
approximations  which  the  results  show  to  be  good  ones.  Finally,  the 
fitted  Gaussian  curves  to  the  entrained  flow  distribution  and  associated 
empirical  formulas  (Table  A10)  are  also  approximations  which  the  results 
show  to  be  good  ones.  These  approximate  formulas  make  the  prediction  of 
entrained  volume  distribution  an  easy  task  within  the  parameter  range 
with  which  this  study  is  concerned. 
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Scales  and  Dimensionless  Ratios  for  Horizontal 
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APPENDIX  B:  COMPARISON  OF  PRESENT  EXPERIMENTS  WITH 
THE  PREDICTIONS  OF  FAN  AND  BROOKS 


Introduction 

1.  By  the  use  of  integral  equation  models.  Fan  and  Brooks  (1969)* 
have  produced  predictions  of  the  behavior  of  buoyant  round  jets  in  a 
linearly  stratified  environment.  Their  model  is  described  in  Appen¬ 
dix  C.  The  relationship  of  these  predictions  to  the  present  work  is 
discussed  in  this  appendix. 

2.  Fan  and  Brooks  expressed  their  results  in  terms  of  different 
parameters  than  those  used  here.  The  conversion  between  their  parameters 
and  the  length  scales  used  in  this  study  is  presented  in  paragraphs  3-13 
of  this  appendix.  A  comparison  of  one  of  the  experiments  with  values 
predicted  by  Fan  and  Brooks  is  given  in  paragraphs  14-15,  and  unified 
diagrams  for  jet  volume  flux  which  incorporate  the  results  of  both  Fan 
and  Brooks  and  the  present  experiments  are  given  in  paragraphs  16-20. 
Concluding  remarks  are  presented  in  paragraphs  21-22. 

Parameter  Conversion 

3.  Fan  and  Brooks  expressed  their  results  as: 


tt,  Ht,  Mt  =  f(m0,  Mc) 


where  »  1t  ,  and  are  terminal  rise  height,  horizontal  distance, 

and  volume  flux  parameters,  and  m^  and  are  source  momentum  and 

volume  flux  parameters  defined  by: 


„  .  (1  i  X-)  s 


2  2 

_  (1  +  AZ)FZ 


(l  «  X2)5 


M°  "  2«,/2\3/2T5/8 


where  X  and  Z  are  the  x  and  z  coordinates  of  the  jet  at  its 
m  t 

terminal  height  (X^  and  Z^  are  equivalent  to  Fan  and  Brooks’  X^  and 
Y^) .  The  term  A  is  a  form  of  the  turbulent  Schmidt  number,  a  constant, 
and  a  is  the  entrainment  coefficient,  also  assumed  constant.  The  term 
is  the  terminal  center-line  dilution  and  6  is  a  scale  factor;  F 
is  a  densimetric  Froude  number: 


and  T  is  a  stratification  parameter: 


where  Ap  =  p 
o  a 


is  the  source  density,  p is  the  ambient 


density  at  the  elevation  of  the  source,  u  is  the  jet  velocity,  and  d 
is  the  nozzle  diameter. 


4.  The  problem  is  to  express  s»o  ,  »  4t  »  0t  »  and  Pt  i 

terms  of  the  variables  and  .  We  begin  by  writing  Equa 


in 


tions  66  and  B7  as: 


m  =  KfV’ 

o  1 


(BIO) 


and 


M  =  K  F1/4,T-5/8 

o  Z 


(Bll) 


where 


*  -  O  *  Q 

R1  4 

1  4A* 


(B12) 


and 


5/8 


(  =  {?.+ 

2  2a1/2x3/2 


(B13) 


5.  Now,  from  Equations  A1  through  A8: 

•G) 


i .  fsVVv*/*  (- 

f  m3/4  V4/  g  "a  V  dz/ 

e  n 


d1/20-1/2 


(B14) 


and 


£  ■  ^  ■  «-/4py4  (-  «»» 


6.  We  now  combine  Equations  BIO  and  Bll  as: 


-3/8 

u  m 

o  o 


*(■ 


k„f1/4t'5/8)(k1f2t"1) 


•3/8 


(B16) 


-3/8p-l/2T-l/4 


and  substitute  Equations  B8  and  B9  in  Equation  B16  to  yield: 


^:3/8  ■  v;3/v'V/4(- 

7.  Similarly,  Equation  BIO  can  be  written: 


*r  - 


1/4 


=  kJ'V'V1'4 


(B17) 


__KuVuVj^_fzy'\m^y,2  (BI8) 

8.  Comparison  of  Equation  B14  with  B17  and  Equation  B15  with  B18 
yields : 


(!) 


1/4  K 


3/8 
1 


-3/8 

u  m 
o  o 


(B19) 


and 


M  _  .,-1/4  1/4 

r  -  K1  m0 

e 


(B20) 


9.  To  complete  the  conversions  between  Fan  and  Brooks'  parameters 
and  our  length  scales,  values  of  Kj  and  are  required.  These  are 

given  in  Fischer  et  al.  (1979),  pp  325  and  371,  as: 


a  -  0.0535  ±  0.0025  for  jets  (B21) 

a  =  0.0833  ±  0.0042  for  plumes  (B22) 

and 


A  =  1.2 


for  jets  and  plumes 


(B23) 


written;  by  coabining  Equations  B30,  B31,  B24,  B25,  A6,  and  A8  we  obtain: 


M 


1/4  -3/8 

3 ?4  =  °' 416  Vo 


for  jets 


(B32) 


and 


ye 

M 


1/4 

3/4 


=  0.518y  m"3/8 
t  o 


for  plumes 


(B33) 


12.  The  thickness  of  the  jet  at  collapse  W  can  be  estimated  by 

la 

assuming 

WL  =  2  wt  (B34) 

where  is  the  terminal  half-width  of  the  jet  defined  by  two  standard 

deviations  from  the  mean  of  the  assumed  Gaussian  distributions.  From 
Fan  and  Brooks'  Equation  73: 

wfc  =  6w^d  =  62  -J2  ay  ■  *^4d  (B35) 

for  horizontal  jets  (when  h  =  mo).  Combining  Equations  B21,  B22,  B34, 
and  B35  yields: 

WT  =  0.3036u  m'1/4d  for  jets  (B36) 

L  to 

and 

W  =  0.47l6y  m  ^4d  for  plumes  (B37) 

la  tO 

13.  The  range  of  parameters  for  which  solutions  are  given  by  Fan 


and  Brooks  is: 


0.001  <  Li  <1 
o 


0.0006  <io  <6 
o 


which  corresponds  (using  Equations  B24-B26)  to: 


1.06  x  10"4  <  ^  <  4.18 

e 


and 


0.21  <  y  <  2.1 

e 


Example 

14.  For  an  example,  we  compare  the  results  of  Experiment  E12  (see 
Tables  A1  and  A2)  with  the  predictions  of  Fan  and  Brooks.  For  this  ex¬ 
periment: 

j9  =  0.118  ^  =  1.87 

z  e 

Using  Equations  B2,  B3,  B24-B30,  B36,  B37,  and  Fan  and  Brooks’  Fig¬ 
ures  14,  15,  and  16,  we  can  evaluate  all  parameters  and  compare  them 
with  the  experimental  results  as  shown  in  Table  Bl. 

15.  As  would  be  expected,  the  results  for  this  high  momentum  jet 
show  better  agreement  with  the  predictions  based  on  jet  coefficients 
than  on  plume  coefficients.  It  would  not  be  expected  that  Fan  and 
Brooks’  predictions  would  be  very  close  to  the  present  results  as  their 
model  was  not  intended  for  application  to  low-buoyancy  jets.  Neverthe¬ 
less,  the  comparisons  shown  in  Table  Bl  suggest  an  encouraging  consis¬ 
tency  between  their  approach  and  ours. 


Total  Volume  Flux 


16.  By  combining  the  predictions  of  Fan  and  Brooks  with  the  present 
experimental  results,  we  can  now  predict  the  volume  flux,  or  dilution,  of 
a  horizontal  buoyant  jet  over  the  whole  parameter  range,  from  momentum- 
dominated  neutrally  buoyant  to  buoyancy-dominated  with  negligible  momen¬ 
tum.  For  linear  stratification  with  no  momentum  flux,  and  0=0  deg, 
Equation  A15  becomes: 


ye 

M 


1/4 
3?4  = 


(B38) 


The  volume  flux  predictions  given  in  Fan  and  Brooks'  Figure  16  were  con¬ 
verted  to  the  form  of  Equation  B38  using  Equations  B24,  B26,  and  B32, 
and  were  combined  with  the  present  results  (Figure  A10)  as  shown  in  Fig¬ 
ure  B1 . 

17.  The  results  clearly  demonstrate  the  asymptotic  behavior  dis¬ 
cussed  in  Appendix  A.  That  is,  as  the  nozzle  size  becomes  small  such 
that  Aq/££  <<  1  •  the  results  become  independent  of  ^q/^e  •  Equa¬ 
tion  B38  then  becomes: 


M 


1/4 

374  = 


(B39) 


The  asymptotic  solution  to  Equation  B39  for  a  buoyancy-dominated  flow, 
i.e.,  a  pure  plume,  can  also  be  simply  decided.  In  this  case: 

M  =  f(B,£>  (BAO) 


which,  following  dimensional  analysis  becomes: 


SB 


3/4t-5/8 


(B41 ) 


where  C 


is  an  experimental  constant.  Equation  B41  can  be  written  in 


the  present  notation  of  Equation  B39  by  combining  it  with  Equations  A7 
and  A8  to  yield: 


M 


1/4 

374 


(B42) 


18.  Although  is  a  fundamental  constant,  we  only  have  Fan  and 

Brooks'  numerical  predictions  from  which  to  estimate  its  value.  Using 
Figure  Bl,  we  would  arrive  at  an  estimate  of  of  1 . 1 ;  thus  Equation 

B42  becomes : 


ye 
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1/4 

3/4 


(B43) 


Figure  Bl  was  derived  using  a  jet  entrainment  coefficient,  however,  and 
a  plume  entrainment  coefficient  is  higher  than  that  of  a  jet.  Therefore, 
the  derivation  of  Figure  Bl  was  repeated,  using  a  plume  entrainment  coef¬ 
ficient,  a  =  0.0833  ,  as  recommended  by  Fischer  et  al.  (1979),  to  yield 
a  value  of  =  1.4  .  Therefore,  Equation  B42  becomes: 


ye 
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1/4 

374 


(B44) 


19.  The  results  can  now  be  plotted  in  the  form  of  Equation  B39  to 

predict  the  behavior  of  horizontal  round  buoyant  jets  of  small  nozzle 

size  over  the  whole  parameter  range,  from  momentum-  to  buoyancy- 

dominated.  This  is  done  in  Figure  B2,  where  the  present  experimental 

results  (Figure  A10,  for  £  /£  <  0.2)  and  the  asymptotic  solutions  for 

v  £ 

a  pure  jet  (Equation  A41)  and  a  pure  plume  (Equation  B44)  are  also  shown. 
Predictions  from  Fan  and  Brooks'  model  (Figure  Bl)  and  the  plume  solution 


(Equation  B43),  both  using  the  jet  entrainment  coefficient,  are  also 
shown.  A  better  approach  to  the  numerical  modeling  would  be  to  use  a 
variable  entrainment  coefficient  dependent  on  local  conditions,  as  recom¬ 
mended  by  List  and  Imberger  (1973).  This  refinement  would  hardly  affect 
the  present  discussion,  however. 

20.  The  predictions  of  Fan  and  Brooks  appear  to  be  consistent  with 
the  present  results.  The  deviation  at  £^/££  2  is  probably  caused 

by  the  breakdown  of  the  assumptions  of  the  Fan  and  Brooks  model  for 
weakly  buoyant  jets.  The  results  of  Figure  B2  suggest  that,  for 
«  and  ,  the  flow  is  plumelike  for  £^/££  <  ~0-8  and  jetlike 

for  £^/££  >  ~2  .  In  between,  both  momentum  and  buoyancy  influence  the 
flow.  If  the  nozzle  is  large  enough  that  £^  is  not  much  smaller  than 
or  ££  ,  then  Figure  B1  should  be  used. 


Concluding  Remarks 


21.  Discussion  of  the  present  problem  in  terms  of  competing  length 
srnles  is  seen  to  be  very  convenient.  In  cases  in  which  the  source 
volume,  buoyancy,  and  momentum  are  all  important,  all  the  length  scales 
expressed  by  Equations  A6,  A7 ,  and  A8  should  be  considered.  The  effect 
of  the  nozzle  geometry  is  confined  to  a  distance  of  several  nozzle 
diameters  from  the  source,  however,  and  is  frequently  of  secondary  impor¬ 
tance  compared  to  the  other  length  scales. 

22.  The  present  experimental  and  theoretical  results  of  Appen¬ 
dixes  A  and  B  show  this  clearly.  If  £^  <<  anc*  >  the  len8th  of 

the  potential  core  of  the  jet  is  much  less  than  the  collapse  length. 

The  jet  can  then  be  analyzed  as  if  generated  by  a  point  source  of  momen¬ 
tum  only.  Even  for  larger  values  of  £  /£  >  the  results  are  not  sensi- 

v  £ 

tive  to  £  /£  ,  with  only  the  volume  flux  parameter  (J  showing  much 

ye 

variation . 
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APPENDIX  C:  INTEGRAL  BUOYANT  JET  MODEL  OF  FAN  AND  BROOKS 


1.  The  integral  model  used  for  strongly  buoyant  jets  in  subroutine 
SBJET  is  based  on  the  model  RBJ  presented  in  Koh  and  Fan  (1970).  The 
material  in  this  appendix  is  extracted  from  their  report.  Note  that  the 
notation  is  slightly  different  from  that  followed  in  the  rest  of  this 
report.  The  model  of  Koh  and  Fan  allows  for  the  merging  of  individual 
round  buoyant  jets  to  a  two-dimensional  flow  similar  to  that  issuing 
from  a  slot  jet.  As  SBJET  may  in  the  future  be  extended  to  apply  to 
merging  jets,  the  round  and  slot  flow  formulations  are  presented  below. 

2.  Consider  a  jet  oriented  at  angle  0q  to  the  horizontal  issuing 
fluid  of  density  p^  and  temperature  T^  into  an  ambient  of  density 
stratification  p  (y)  and  temperature  stratification  T  (y)  .  Let  Q 

be  the  discharge,  the  momentum  flux,  Fj  the  density  deficiency 

flux,  and  the  temperature  deficiency  flux  at  the  source.  Figure  Cl 

illustrates  the  general  behavior  of  such  a  jet.  The  bending  of  the  jet 
path  is  a  result  of  the  fact  that  the  discharge  is  buoyant.  Define  u* 
as  the  velocity,  T*  as  the  temperature,  and  p*  as  the  density  of  the 
buoyant  jet  flow.  Since  the  ambient  is  motionless,  u*  is  assumed  to 
be  along  the  jet  path.  Let  s  be  the  coordinate  along  the  jet  path, 
A-plane  be  the  plane  perpendicular  to  the  jet  path,  and  0  the  angle 
of  the  jet  path  with  respect  to  horizontal.  We  now  define  the  volume 
flux  Q  ,  momentum  flux  M  ,  density  deficiency  flux  F  ,  and  tempera¬ 
ture  deficiency  flux  G  along  the  jet  trajectory  as: 


Q  =  J u*dA 
A 


i  f  *2  *  r 

—  jf  u  p  dA 


o  A 


(pa  -  p*)u*dA 


*2 

u  dA 


(Cl) 

(C2) 

(C3) 


ro 


Figure  Cl.  Definition  diagram  for  integral  model 
of  round  buoyant  jets 


G  = 


T*)  u*dA 


(CA) 


Note  that  M'  is  the  true  momentum  flux  while  M  is  the  kinematic 
momentum  flux. 

3.  We  note  that  in  the  vertical  direction,  a  buoyancy  force  exists 
due  to  the  density  difference  between  the  jet  fluid  and  the  ambient  fluid 
tending  to  bend  the  jet.  This  force  f  is 

f  =  g  f  (Pa  -  P*)dA  (C5) 

A 

A.  The  conservation  equations  can  now  be  written  in  terms  of 
these  variables.  The  conservation  of  mass  equation  is 
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C2 


AQ  =  EAS 

where  E  is  the  rate  of  entrainment  of  ambient  fluid.  Note  that, 
strictly  speaking,  since  the  density  is  variable,  we  should  really  have, 
instead  of  Equation  C6 


[/uVd*]  = 


However,  all  density  differences  are  small,  and  we  may  approximate  Equa¬ 
tion  C7  by  C6.  The  variable  AQ  is  equivalent  to  Q  (see  Equation  11, 
main  text). 

5.  The  conservation  of  horizontal  momentum  flux  is 


d(M  cos  6)  _  Q 
ds 


For  the  vertical  momentum,  we  must  include  the  buoyancy  force.  Thus, 


-5—  (M'  sin  0)  =  f 
ds 


The  conservation  of  density  deficiency  flux  equation  reads 


[{■* 


h  f  U'V(Po  ‘  pk)dA  =  E(Po  •  Pa) 


(CIO) 


where  Pq  is  a  reference  density  (e.g.,  pQ  =  Pa  (0)).  Equation  CIO 
can  be  written 

h  [/■**.  '  Pa)flA  +  /u*(p  a  "  P")dAj  =  E(Po  "  Pa} 
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Using  Equation  C6  we  finally  have 


dF  dPa 


ds  ds 


(C12) 


Similarly,  the  conservation  of  temperature  deficiency  flux  equation  reads 


£  f  “X  -  T*>dAl 


(T  -  T  )dA  =  E(T  -  T  ) 
o  o  a 


CC13) 


which  reduces,  with  Equation  C6,  to 


dG  dTa 


ds  ds 


(C14) 


6.  Equations  C6,  C8 ,  C9 ,  C12,  and  C14  constitute  five  equations 
for  the  five  unknowns  Q  ,  M  ,  0  ,  F  ,  G  ,  as  functions  of  s  , 
once  we  can  express  E  and  f  in  terms  of  known  quantities  or  thpse  un¬ 
knowns.  To  do  this,  we  will  make  two  more  assumptions.  First,  we  shall 
assume  similarity  of  the  shapes  of  the  velocity  profile,  temperature  de¬ 
ficiency  profile,  and  density  deficiency  profile  in  the  plane  A  .  In 
particular,  it  will  be  assumed  that  the  profiles  are  Gaussian.  Thus,  in 
the  two-dimensional  case  (slot  jet)  we  assume 


*  ,  v  -n2/b2(s) 

u  (s,n)  =  u(s)e 


(CIS) 


Pa  -  p*(s,n)  =  [pa  -  p(s)]  e 


-n2M2b2 


wmm 


f  *.  «.  • 


jv  -v  v-v-v 


»•  v  o  vivC 


j-./.-.-.-.v-- 

►  .  .\s  .  . 

.v 


‘T.v7-/v'.y. 


(C16) 


where  u(s)  ,  p(s)  ,  and  T(s)  are  the  values  along  the  jet  center 
line.  The  term  r|  is  the  coordinate  normal  to  s  ,  b(s)  is  the  char¬ 
acteristic  jet  width,  and  X  is  a  turbulent  Schmidt  number  for  the  two- 

s 

dimensional  case  (2-D).  Similarly,  in  the  axisymmetric  case,  we  take 


Pa  -  P 


u  (s  ,  r)  =  u(s)e 
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where  r  is  the  radius  normal  to  the  s  coordinate. 

7.  Secondly,  we  shall  assume  that  the  entrainment  function  E  is 
proportional  to  the  jet  characteristic  velocity  u  and  the  jet  boundary 
(2nb  or  2L  for  the  axisymmetric  and  2-D  case,  respectively)  and  that 
the  proportionality  constant  is  a  (a^  for  round  jet  and  for  slot 

jet).  Substituting  these  expressions  into  the  definitions  for  Q  ,  M  , 
F  ,  G  ,  and  f  (Equations  C1-C5)  will  give  Q  ,  M  ,  F  ,  G  ,  and  f 
expressed  in  terms  of  the  quantities  u  ,  p  ,  T  ,  and  b  .  For  ex¬ 
ample,  substituting  Equation  C15  into  Cl  gives 


oo  2  2 

Q  =  L  /  u(s)e  ^  ^  dij 


where  L  is  the  length  of  the  2-D  slot  jet.  Thus, 
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9.  Table  Cl  gives  the  transformation  from  the  variables  u  ,  b 
etc.,  to  the  variables  Q  ,  M  ,  etc.  (The  inverse  transformation  is 
given  in  Table  C2.  Moreover,  it  is  possible  to  express  E  and  f  now 
in  terms  of  Q  ,  M  ,  etc.,  as  shown  in  Table  C3.) 

10.  The  problem  of  the  mixing  processes  involved  for  a  row  of 
round  buoyant  jets  spaced  a  distance  L  apart  is  now  considered.  Ini¬ 
tially,  the  jets  are  separate  round  jets.  However,  after  a  while,  they 
begin  to  merge  and  form  more  nearly  a  two-dimensional  slot  jet.  Thus, 
in  the  calculations,  it  is  necessary  to  provide  a  criterion  whereby  the 
round  jet  analysis  is  switched  to  that  for  a  slot  jet.  Two  such  cri¬ 
teria  are  proposed.  We  assume  that  transition  occurs  when  the  width  of 
the  round  jet  becomes  equal  to  the  jet  spacing.  This  shall  be  desig¬ 
nated  transition  1.  Referring  to  Table  C2,  this  occurs  when 


or  — =  L  £2  = 
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where  the  "jet  width"  is  taken  as  2  ^ 2b  . 


Table  C2 

Inverse  Transformation  of  Variables,  Buoyant  Jet  Model 


dx  = 

ds 


sin  6 


(C23) 


The  system  of  Equations  C 6,  C8,  C 9,  C12,  C14,  C22,  and  C23  constitute 
seven  ordinary  differential  equations  for  the  seven  unknowns,  Q  ,  M  , 

6  ,  F  ,  G  ,  x  ,  and  y  as  a  function  of  s  .  These  equations  may 
be  solved  given  the  initial  values  of  the  unknowns  at  s  =  0  . 

12.  The  initial  conditions  are  given  by  the  source  conditions, 
namely,  uq  ,  the  jet  velocity;  Dq  ,  the  jet  diameter;  T^  ,  the  jet 
temperature;  ,  the  jet  density;  and  0q  ,  the  jet  discharge  angle. 

However,  since  the  formulation  is  in  terms  of  the  flux  quantities  Q  , 

M  ,  F  ,  and  G  ,  these  jet  characteristics  must  be  converted  to  initial 
values  in  these  variables.  Moreover,  it  is  well  known  that  there  exists 
a  zone  of  flow  establishment  extending  a  few  jet  diameters  during  which 
the  top-hat  profiles  of  velocity,  density  deficiency,  and  temperature  ex¬ 
cess  change  gradually  to  Gaussian  form  (Figure  C2).  In  this  formulation, 
we  shall  start  the  integration  from  the  beginning  of  the  zone  of  estab¬ 
lished  flow.  Thus,  it  is  necessary  to  relate  the  jet  characteristics  to 
the  flux  quantities  at  this  point.  Albertson  et  al.  (1950),  in  their  ex¬ 
perimental  investigations  on  the  round  jet,  found  that  the  zone  of  flow 
establishment  extends  a  distance  of  6.2  jet  diameters.  Equating  the  mo¬ 
mentum  flux  at  the  beginning  and  end  of  the  zone  of  flow  establishment, 
assuming  that  the  buoyancy  force  is  negligible  in  such  a  short  region,  we 
get  (see  Table  Cl  for  M  and  round  jet) 
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where  b^  is  the  jet  half-width  and  u  equals 
zone  of  flow  establishment. 


u  at  the  end  of  the 
o 


13.  Thus,  the  initial  value  for  Q  is 


In  other  words,  the  volume  flux  at  the  beginning  of  the  zone  of  estab¬ 
lished  flow  is  twice  that  at  the  source. 

14.  By  assuming  further  that  the  ambient  density  is  uniform  in  the 
zone  of  flow  establishment,  we  may  equate  the  density  deficiency  flux  at 
the  beginning  and  end  of  this  zone  to  get 
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Thus,  the  initial  value  for  F  is 


F  =  7  D2u  (p  -  p,) 
1  4  o  o  a  1 


However,  the  center-line  density  deficiency  is 
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Similarly, 
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and  the  center-line  temperature  excess  is 
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15.  The  equations  presented  above  are  solved  numerically  by  the 
program  RBJ  (subroutine  SBJET).  The  method  of  solution  is  to  first  ob¬ 
tain  the  initial  conditions  (Q^ ,F^ ,G^)  from  the  given  source  charac¬ 
teristics.  Then  the  Equations  C6,  C8,  C9 ,  C12,  C14,  C22,  and  C23  are 
integrated  with  E  and  f  as  given  by  those  for  the  round  jet  (col¬ 
umn  2,  Table  C3).  When  transition  is  reached,  given  by  Equation  C21, 
one  simply  continues  the  solution  but  with  E  and  f  as  given  by  those 
for  the  slot  jet  (column  3,  Table  C3).  The  results  obtained  are  then 
converted  from  the  variables  Q  ,  M  ,  F  ,  G  ,  etc.,  to  the  physical 
variables  u  ,  p  ,  T  ,  and  W  ,  the  jet  width  that  is  taken  to  be 

2  2b.  The  conversion  is  effected  by  the  relations  in  Table  C2. 

16.  the  input  to  the  program  consists  of  the  following: 

u  =  jet  velocity 
o 

D  =  jet  diameter 
o  J 

Tj  =  jet  temperature 


Cll 


m  % 


t\ 


jet  density,  gm/cm 

jet  discharge  angle  with  respect  to  horizontal,  deg 
jet  discharge  depth 
jet  spacing 

entrainment  coefficient  for  a  round  jet 
entrainment  coefficient  for  a  slot  jet 
spreading  ratio  for  a  round  jet 
spreading  ratio  for  a  slot  jet 
gravitational  acceleration 


